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INTRODUCTION 


QuaHty  control  (QC)  should  be  an  integral  part  of  any  data  coUection 
program.  The  impetus  to  perform  quality  control  is  derived  from  management’s 
quality  assurance  (QA)  program.  The  QA  program  is  the  collection  of  plaimed 
and  systematic  actions  defined  by  management  as  those  necessary  to  provide 
confidence  that  a  product  or  service  will  satisfy  customer  quahty  requirements. 
This  program  should  be  defined  in  management’s  quality  policy  and  imple¬ 
mented  through  a  quality  system  that  includes  the  organizational  structure, 
responsibilities,  procedures,  processes,  and  resources  for  implementing  qu^ity 
management.  Quality  control  includes  the  operational  techniques  and  activities 
needed  to  fulfill  QA  program  requirements  (ANSI/ASQC,  1987).  Further 
impetus  towards  quality  assimance  comes  from  increased  scrutiny  of  and 
accountability  for  environmental  impacts  of  range  activities  and  the  need  to 
meet  International  Standards  Organization  (ISO)  standards.  This  guide  is 
intended  to  assist  ranges  in  the  development  and  maintenance  of  an  acceptable 
quality  control  program  for  surface  meteorological  data. 


IX 


CHAPTER  1 


QUALITY  CONTROL  GENERAL  PRINCIPLES 


1.1  Product  and  QuaRtv  Control 

The  range  meteorologist’s  job  is  to  provide  valid  and  relevant  information  to 
the  customer.  To  achieve  this  goal,  the  information  must  be  supportable  by 
documentation.  Validity  can  be  defined  using  instrument  performance  and 
calibration  records  and  a  description  of  the  steps  used  to  process  the  data.  Rele¬ 
vance  which  includes  the  timely  dissemination  of  representative  and  pertinent 
infoniation  in  a  user-friendly  format,  is  as  important  as  validity  to  a  successful 
measurement  program.  Collecting  valid  but  irrelevant  data  contributes  nothmg  to 
program  objectives  and  wastes  resources.  To  provide  the  customer  vdth  a  quality 
product,  a  measurement  program’s  quality  control  (QC)  plan  should  include 
procedures  designed  to  meet  the  program’s  requirements  in  terms  of  validity  and 
relevance  and  to  present  this  valid,  relevant  information  in  an  accessible  and 
usefiil  format.  Consequently,  QC  must  begin  with  a  dear  definition  of  program 
objectives  that  can  be  satisfied  by  a  specific  product  set.  Product  definition  is 
followed  by  careful  identification  of  relevant  measurement  requirements, 
instrumentation,  and  sites. 

An  optimiun  QC  program  indudes  proactive  identification  and  avoidance  of 
potential  errors  as  well  as  error  checking  procedures.  The  QC  process  starts 
before  a  measurement  program  and  might  proceed  as  described  in  the  following 
subparagraphs. 

1.1.1  Reauirementfi  DpAnition.  Detailed  information  regarding  the  measurements 
(and  allowable  measurement  uncertainties)  needed  to  provide  that  information, 
and  the  instrumentation  used  to  obtain  those  measurements. 

1.1.2  Calibration/Certification.  A  measurement  comparison  made  against  ^ 
accepted  standard  or  transfer  standard  or  against  some  standard  in  a  functional 
test,  as  appropriate. 

1.1.3  Tngfmmpnt  and  Site  Selection.  Choosing  instruments  and  measurement 
sites  that  provide  representative  measxirements  of  pertinent  variables. 

1.1.4  Meteorologif»a1  T.imitations.  Go/no-go  controls  or  restrictions  imposed  upon 
the  mission  or  measmement  program. 
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1.1.5  Deployment.  Installation  of  the  instruments  and  completion  of  field  checks 
and  reasonableness  (consistency)  checks. 

1.1.6  Flagging.  Identification  of  out-of-range  or  out-of-tolerance  conditions. 

1.1.7  Monitoring.  Real-time  review  of  incoming  data  as  needed  to  satisfy  mission 

requirements. 

1.1.8  Analysis.  Review  of  summary  statistics  and  plots  derived  from  the  data  set. 
Performing  flagging/de-flagging,  spike  removal,  and  trend  analysis. 

1.1.9  Report.  Preparation.  Conversion  of  the  data  into  meaningful  information 
that  is  supported  by  documentation  of  the  instrumentation,  measurement,  and 
analysis  processes. 

1.1.10  ModeUng.  Applying  the  information  through  prototypes  to  generate  new 
knowledge  or  understanding  of  meteorological  processes. 

The  goal  of  QC  is  to  provide  the  desired  product  to  both  the  extern^  customer  who 
funds  the  program  and  to  the  internal  customer  who  is  the  performing 
organization  itself. 

1.2  Experimental  and  Mission  Support  Process^ 

The  activity  of  generating  meteorological  information  generally  follows  one 
of  two  general  processes:  (1)  the  Experimental  Process,  or  (2)  the  Mission  Support 
Process.  Each  process  consists  of  a  series  of  actions  and  products.  Quality 
Control  is  embedded  within  each  of  the  actions  taken  to  produce  the  desired  prod¬ 
ucts.  Some  QC  procedxu-es  are  more  appropriate  for  one  process  &an  for  the 
other.  Selecting  the  appropriate  QC  procedures  for  a  process  is  vital  to  the  design 

of  a  QC  program. 

The  Experimental  Process  begins  with  a  question  that  leads  to  a  test  or 
experiment  that  produces  an  initial  product  (data).  The  data  are  vahdated  and 
analyzed  to  produce  information  as  the  next  product.  The  process  continues 
through  modeling  to  produce  knowledge,  which  serves  as  a  basis  for  hypothesizing 
new  questions  or  theories  that  lead  to  further  testing  and  the  generation  of  new 
data  sets.  The  continuous  aspect  of  this  process  has  QC  inherently  embedded 
within  it;  errors  surface  as  inconsistencies  that  trigger  review,  re-an^ysis, 
corrective  action,  and  re-testing.  The  Experimental  Process  can  be  viewed  as  an 
ongoing  ascending  spiral  of  actions  and  products.  The  upper  half  of  figure  1-1  is  a 
two-dimensional  representation  of  the  Experimental  Process  spiral  with  actions 
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Figiire  1-1.  The  Experimental  Process  and  the  Mission  Support  Process. 


illustrated  as  arrows  and  products  illustrated  as  boxes.  This  figure  describes  &e 
process  of  scientific  inquiry  as  it  might  be  performed  in  university  and  research 
settings.  Data  collection,  analysis,  and  archival  for  climatological  purposes  may 
also  follow  the  Experimental  Process. 

Meteorological  data  obtained  in  the  Experimental  Process  are  as  inclusive 
as  possible  within  experimental  and  budgetary  limitations.  The  majority  of  the 
QC  is  performed  well  after  data  coUection  has  been  completed  when  there  is  time 
for  reflective  analysis  and  inteipretation.  Data  are  carefully  formatted,  annotated, 
and  archived  for  future  use.  Unexpected  results  such  as  anomalous  reading  or 
outlier  data  points  are  often  the  subject  of  intense  scrutiny.  Quahty-control 
procedures,  which  may  include  the  use  of  complex  statistical  or  analytical 
techniques,  are  reported  in  detail  for  peer  review,  where  the  data  and  their 
processing  procedures  are  subject  to  re-analysis  and  alternate  interpret^ons. 

The  Experimental  Process  employs  the  principles  of  the  scientific  method  to 
identify  and  to  preclude  the  retention  of  error.  The  discovery  of  error  often  leads 
to  new  measure-magnet  methods  or  experimental  procedures  designed  to  elimmate 

those  errors. 

The  Mission  Support  Process  more  closely  resembles  the  activities  at  test 
ranges  where  real-time  go/no-go  decisions  must  be  made  within  a  time  constraint. 
Unlike  the  Experimental  Process,  the  Mission  Support  Pro^ss  is  driven  by  user- 
defined  start  and  end  times.  These  times  are  often  determined  by  factors  beyond 
the  control  of  the  data  analyst,  who  must  produce  timely  information  on  which 
critical  mission  decisions  will  be  based.  The  lower  half  of  figure  1-1  illusteates  the 
Mission  Support  Process.  Data  collected  for  mission  support  ^e  analyzed  and 
validated  to  the  point  where  information  is  available  for  decision  making  and  not 
necessarily  to  the  point  where  sufficient  knowledge  accimulates  to  identify  all 
errors  or  to  take  cogent  corrective  actions.  The  emphasis  is  on  producing 
consistent  information  to  support  (or  refute)  existing  expectations  (forecasts)  of 
meteorological  conditions.  Outhers  must  be  either  accepted  as  unanticipated 
results  or  rejected  as  erroneous  data  based  on  the  judgement  of  the  on-duty 
analyst.  Measurement  redundancy  is  often  needed  to  ensure  a  reasonable  degree 
of  success  in  having  the  required  valid  information  available  before  decision  dead¬ 
lines.  Data  validation  often  is  based  on  rule-based  "eiffier/or"  criteria  rather  tl^ 
detailed  statistical  analyses.  Mission-support  information  is  usually  presented  in 
brief  statistical  summaries,  while  the  raw  data  are  stored  in  unique  formats  that 
are  difficult  to  use  outside  the  originating  test  center.  Archival  is  not  performed 
with  a  view  to  re-analysis  or  for  peer  review,  so  the  data  often  do  not  resurface 
unless  something  goes  wrong  with  the  mission.  The  Mission  Support  Process, 
unlike  the  Experimental  Process,  is  not  self-checking.  A  conscious  effort  is  needed 
to  build  quality  checks  into  the  data  collection  and  analysis  phases  of  the  Mission 
Support  Process.  This  guide  is  designed  to  assist  the  range  meteorologist  with 

that  requirement. 
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1.3  Validation 


Meteorological  measiirement-system  validation  includes  some  combination 
of  calibration,  on-site  system  checks,  and  functional  testing.  Its  purpose  is  to 
assme  the  user  that  the  measurement  system  is  performing  to  a  known  standard. 
Calibrations  usually  verify  system-component  response  over  a  desired  range  or 
span,  while  system  checks  verify  the  integrity  of  signal  transfer  through  the 
measurement  system.  Functional  testing  includes  a  statistical  method  to  compare 
the  performance  of  instruments  measuring  the  same  phenomena  xmder  field 
conditions.  Results  of  calibrations,  system  checks,  and  functional  testing  are 
expressed  using  terminology  presented  in  table  1-1.  Documentation  (log  records) 
of  calibration,  maintenance,  and  system  checks  is  needed  to  defend  measurement 

validity. 

13  1  Validation  and  Accuracy.  System  validation  usually  includes  a  statement  of 
performance  in  terms  of  comparability  or  accuracy.  This  statement  should  address 
the  entire  system  not  just  component  parts.  Total  measurement-system  accuracy 
is  the  product  of  uncertainties  attributable  to  the  sensor,  electronic-data  tra^- 
mission,  and  data-logging  process.  Sensor  response  includes  an  analog  or  digit^ 
signal  derived  from  the  sensor’s  interaction  with  the  environment.  This  signal  is 
then  subject  to  line  losses  or  noise  during  transmission  to  a  recording  device.  The 
recorder  also  contributes  its  noise  threshold  to  total  system  noise.  Measurement 
accuracy  is  maximized  by  choosing  a  compatible  sensor/recorder  combination, 
careful  grounding  and  cabling  procedures,  filtering,  and  flagging  suspect^  data  for 
review.  Independent  errors  occurring  at  each  step  in  the  process  result  in  a  total 
error  that  is  equal  to  the  square  root  of  the  sum  of  the  squares  of  each 
contributing  error. 

1.3.2  rifllibrfltinn.  Instrument  calibration  is  the  first  step  in  defining  data 
validity.  Calibration  involves  comparison  against  a  known  standard  to  determine 
how  closely  instrument  output  matches  the  standard.  Performance  of  a  laboratory 
calibration  carries  the  imphcit  assumption  that  the  instrument  s  diaracteristics 
are  sufiSdently  stable  to  retain  that  calibration  in  the  field.  A  ^bration  history 
established  over  successive  calibrations  can  provide  confidence  in  the  instrument’s 
stability. 

Laboratory  calibration  may  define  a  bias  which  can  then  be  removed 
through  mechanical,  electrical,  or  software  adjustment.  The  remaining  random 
error  or  imprecision  is  not  repeatable  and  cannot  be  removed,  but  it  can  be 
statistically  defined  through  a  sufficient  number  of  measurement  repetitions 
during  cahbration.  The  combined  effect  of  all  random  (uncorrelated)  uncertainties 
can  be  obtained  from  the  square  root  of  the  sum  of  the  squares  of  the  individual 
uncertainties. 
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TABLE  1-1.  DEFINITION  OF  TERMS 


,  . . . 

Accuracy:  The  absolute  imcertainty  in  a  measurement  obtained  by  comparison 
against  an  accepted  standard  in  a  carefully  controlled  environment;  the  degree 
to  which  a  measurement  corresponds  to  a  "true"  value  which  is  never  absolutely 
known.  In  practical  applications,  an  instrument’s  accuracy  is  determined  by 
comparison  of  its  measurements  to  corresponding  measurements  made  by  a 
"standard"  instrument  with  differences  expressed  as  a  "bias"  plus  or  minus  a 
random  uncertainty  or  precision. 

Bias:  The  mean  difference  between  averaged  readings  of  the  test  instrument 
and  the  standard  instrument  to  which  it  is  compared. 

Comparability:  The  root  mean  square  (rms)  of  the  difference  between  readmgs 
of  two  or  more  instruments  of  different  design  sampling  the  same  variable  in  an 
uncontrolled  field  environment. 

Functional  Precision:  The  rms  of  the  difference  between  readings  of  two  or 
more  instruments  of  the  same  design  sampling  the  same  variable  in  the  same 
environment. 

Precision:  Response  of  the  measuring  instrument  to  changes  in  the  measured 
variable.  Components  of  precision  include  resolution  and  repeatability. 

Repeatability:  The  ability  of  an  instrument  to  return  to  a  known  measurement 
or  position  after  being  cycled  through  its  operating  range. 

"Rjinprp-  The  span  of  measurement  which  extends  fi:om  the  threshold  to  the 
mflYimiiTn  measurable  value. 

Resolution:  The  least  significant  digit  reading  which  defines  the  smallest 
variation  in  the  vEiriable  input  that  causes  a  detectable  change  in  instrument 
output.  It  is  given  as  a  fi*actional  part  of  the  measurement  scale. 

ffamnlinfr  Interval:  The  period  between  successive  instrument  readings. 
Considerations  are  given  to  sample-size  requirements,  instrument-time 
constant,  data-logging  speed  and  capacity,  and  measurement  influence  sampling 
interval  purpose. 

Threshold:  The  smallest  measurable  input. 

Timp  Cnnstflnt-  A  measure  of  an  instrument’s  rate  of  response  to  a  step  change 
in  the  measured  variable.  It  is  the  time  required  for  the  instrument  response  to 
reach  |l-l/e|,  or  63  percent,  of  the  new  equilibrium  value. _ 


Calibration  should  define  an  instrument's  bias  or  average  deviation  fi-om  its 
reference  standard,  the  range  over  which  the  calibration  is  valid,  and  the  existence 
of  any  thresholds  or  nonlinear  response  regions.  It  should  also  define  an 
instrument’s  resolution,  which  is  ^e  smallest  change  in  input  that  produces  a 
measurable  change  in  output.  Hysteresis,  which  represents  the  instrument  s 
imperfect  abihty  to  gain  or  shed  energy,  should  be  identified  by  cychng  the  sensor 
over  its  operating  range  during  calibration. 

The  pedigree  of  a  calibration  depends  on  its  traceability  (MIL-STD-45662A, 
1980).  Traceability  is  defined  by  a  hierarchy  of  standards.  The  primapr  standard 
is  the  ultimate  reference  with  respect  to  fundamental  physic^  quantities  (mass, 
length,  time)  and  offers  the  highest  obtainable  precision.  Primary  standards 
reside  within  major  national  or  international  institutions.  Primary  standard 
instruments  are  usually  impractical  for  field  measurements.  Secondary  standards 
are  traceable  to  a  primary  standard.  These  instruments  often  reside  in  major 
calibration  laboratories  and  are  usually  impractical  for  field  use.  Working 
standards  are  usually  laboratory  instruments  that  have  been  calibrated  against  a 
secondary  standard.  Working-standard  instruments  that  are  actually  taken  to 
nonlaboratoiy  sites  are  known  as  transfer  standards.  Transfer  standard  instru¬ 
ments  be  used  to  compare  instruments  in  a  laboratory  (Lockhart,  1987)  or  the 
field  (Lockhart,  1989B). 

Traceability  to  a  recognized  standard  is  available  for  some  meteorological 
instruments.  The  National  Institute  for  Science  and  Technology  (NIST)  (formerly 
National  Bureau  of  Standards)  maintains  a  piston-gage  standard  for  pressure,  an 
International  Practical  Temperature  Scale  (IPTS-90)  for  temperature,  low-  and 
high-speed  wind  tunnels  for  wind  speed,  and  a  gravimetric  device  for  humidity 
(Brock  and  Nicolaidis,  1984).  The  World  Meteorological  Organization  (WMO) 
mainfains  an  absolute  cavity  pyrheliometer  at  Davos,  Switzerland,  with  a 
secondary  National  Oceanic  and  Atmospheric  Administration  (NOAA)  standard  in 
Boulder,  Colorado.  When  no  recognized  standard  is  available,  a  consensus 
standard  may  be  established  to  provide  a  measure  of  relative  agreement  between 
instruments  (MIL-STD-45662A). 

When  possible,  it  is  desirable  to  have  a  calibration  referenced  by  secondary 
or  transfer  standards  to  NIST  and  measurement  procedures  referenced  to  a 
standards  organization  such  as  the  American  Society  for  Testing  and  Materials 
(ASTM)  that  is  approved  by  the  American  National  Standards  Institute  (ANSI). 
These  standards  are  available  for  some  types  of  instrumentation,  but  for  others, 
especially  those  representing  newer  technologies,  they  are  largely  unavailable. 
Standards  are  most  often  available  for  instruments  that  produce  analog  voltage  or 
current  outputs  that  can  be  checked  statically  from  mounting  jigs  or  isothermal 
baths  or  dynamically  in  a  wind  tunnel.  In  general,  industry  standards  have  not 
kept  up  with  advances  in  remote  sensing  and  microprocessor  technology. 
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Microprocessors  typically  perform  some  degree  of  processing  or  mterpreta- 
tion  of  a  sensor’s  output,  which  complicates  calibration.  In  the  absence  of  iMustry 
standards,  an  instrument’s  manufacturer  should  provide  cahbration  and  perfor¬ 
mance  guidance.  Initial  cahbration  cycles  may  require  adjustment  because  of  the 
instrument’s  exposure  to  harsh,  corrosive  environments  or  severe  weather  events. 

A  logbook  or  record  of  an  instrument’s  cahbration  history,  supported  by 
whatever  cahbration  documentation  is  available,  should  be  kept  on  file  for 
reference  and  should  accompany  each  instrument  returned  to  a  laboratory  for 
cahbration.  Each  cahbration  log  should  include  the  instrument’s  departure  fi:om 
its  standard  as  it  enters  the  laboratory  for  cahbration.  The  cahbration  deviation 
can  then  be  compared  with  the  instrument’s  measurement  accuracy  requiremeiits. 
Cahbration  intervals  should  initiaUy  follow  manufacturer’s  reco^endations.  An 
instrument  that  establishes  a  record  of  deviations  that  fah  within  acceptable 
accuracy  limits  may  be  returned  for  calibration  at  less  firequent  intervals. 
Conversely,  an  instrument  that  consistently  arrives  at  a  cahbration  facility  in  an 
out-of-tolerance  condition  may  require  replacement  or  more  fi*equent  cahbrations. 


The  NIST  traceabihty  cannot  be  the  sole  basis  for  evalimting  the  validity  of 
an  instrument’s  reading.  Instrumentation  for  which  traceabihty  to  a  NIST 
primary  or  secondary  standard  has  not  been  established  can  produce  vahd  data  if 
referenced  to  a  fimdamental  principle.  For  example,  no  NIST  standard  exists  for 
a  sonic  anemometer/thermometer,  but  speed-of-sound  propagation  is  fimdamen- 
tally  related  to  the  temperature  of  the  propagation  medium.  The  propagation  of 
an  acoustic  wavefiront  in  a  stih  (zero  wind)  chamber  contmning  a  pure  gas  at  a 
known  temperature  provides  a  fundamental  basis  for  cahbration.  A  zero-wind 
check  in  air  provides  a  method  to  verify  calibration  stabihty  in  the  field. 
Conversely,  a  manufacturer’s  claim  that  a  sensor  is  NBS  or  NIST  traceable  may 
be  of  httle  value  if  it  leads  to  an  accepted  standard  through  a  long,  undocumented 
cahbration  trail  with  compounded  imcertainties.  The  NIST  cannot  guarantee  the 
accuracy  or  precision  of  a  commercial  instrument  even  if  it  is  described  as  NIST 
traceable."  Traceabihty  can  be  estabhshed  through  an  accredited  calibration 
laboratory  such  as  the  Test  Measurement  and  Diagnostic  Equipment  (TMDE) 
Facility  at  Redstone  Arsenal. 

13  3  Fimctional  Testing.  Calibration  provides  a  necessary  but  not  sufficient 
basis  for  defining  data  accuracy,  because  cahbration  of  an  instrument  in  a 
laboratory  cannot  completely  define  how  weU  it  wiU  operate  in  the  field.  For 
example,  a  well-calibrated  propeller  anemometer  operated  in  its  stah  zone  (wind 
flow  within  a  few  degrees  of  perpendicular  to  the  propeller  shaft)  wiU  not  produce 
usable  data.  Fimctional  testing  is  a  method  for  defining  how  well  instruments 
perform  in  the  field,  and  offers  the  advantage  of  comparing  instruments  while 
they  are  performing  imder  real  atmospheric  conditions.  Figures  of  merit  used  in 
functional  testing  are  bias,  comparability,  and  precision  as  defined  by  Hoehne 
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(1971)  and  presented  in  table  1-1.  These  figures  of  merit  are  used  to  define 
performance  with  respect  to  an  accepted  standard  instrument  exposed  to  thy  ame 
conditions.  Functional  precision,  as  defined  by  Hoehne  (1977),  is  the  rms  differ- 
ence  between  readings  of  two  or  more  instruments  of  the  same  design  operating  in 
the  game  environment.  Functional  precision  defines  whether  differences  in 
readings  fi'om  two  similar  instruments  are  likely  to  be  significant.  Differences  m 
instrument  readings  less  than  the  fimctional  precision  are  unlikely  to  be  signif¬ 
icant.  Comparability  is  an  analogous  term  used  when  testing  instruments  of 

different  design. 

For  practical  applications,  good  operational  comparability  or  functional 
precision  between  instruments  is  a  more  valuable  attribute  than  precise  absolute 
calibration.  Functional  testing  can  be  performed  in  the  field  using  a  transfer 
standard  (Lockhart,  1989B),  saving  instrument  downtime.  For  many  meteorolog¬ 
ical  quantities,  accurate  differences  in  instrument  readings  such  as  wind  she^s 
and  temperature  gradients  are  more  important  than  absolute  accuracy  of  individ¬ 
ual  readings.  A  disadvantage  of  fimctional  testing  is  that  the  tested  instruments 
are  rarely  exercised  over  their  full  operational  range.  Functional  testing, 
therefore,  compliments  but  does  not  replace  laboratory  cafibration.  On  the  other 
hand,  fimctional  testing  may  be  the  only  practical  means  available  to  determine 
the  performance  of  remote-sensing  instruments  that  cannot  be  operated  in  a 
laboratory.  Finkelstein  et  al.  (1986)  and  Hoehne  (1977)  provide  examples  of 
fimctional  testing  for  in-situ  meteorological  instruments. 

1.3.4  ffipnfll  ObRcks.  Installation  of  calibrated  analog  signal-generating  eqmp- 
ment  is  followed  by  end-to-end  signal  checks  prior  to  test  data  collection.  Signal 
checks  involve  application  of  a  known  voltage  or  firequency  at  the  sensor  head  and 
monitoring  the  result  at  the  point  where  data  are  logged.  Translators  often 
include  a  voltage-bias  adjustment  designed  to  compensate  for  system  bias.  Alter¬ 
natively,  the  signal  can  be  logged  with  its  known  bias  and  subjected  to  later  soft¬ 
ware  corrections.  Slow  drift  in  signal  bias  because  of  grounding  changes  can  be 
identified  by  periodically  repeating  end-to-end  checks  or  performing  these  checks 
when  significant  changes  in  the  electromagnetic  environment  are  suspected. 
Instruments  that  provide  a  shielded  digital  output  are  less  susceptible  to  signal 
corruption.  Interference  large  enough  to  affect  the  signal  usually  causes  data  loss 

1.4  Relevance 

Relevance  is  the  first  and  final  consideration  in  a  QC  program.  Unlike 
validity,  relevance  is  often  a  qualitative  rather  than  a  quantitative  matter  teat 
reqmres  application  of  judgment  and  communication  skills  as  well  as  tec^cal 
knowledge.  Establishment  of  a  relevant  measurement  program  begins  with  an 
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imderstanding  of  program  objectives  followed  by  selection  of  adequately 
instrument  sites,  the  use  of  appropriate  sampling  and  averaging  times,  and  the 
timely  dissemination  of  the  required  information  in  a  usable  format. 

141  Measurement  Program  Planning.  Relevance  is  particularly  difficult  to 
achieve  for  the  Mission  Support  Process,  because  data  must  be  qmckly  gathered 
and  transformed  into  information  without  the  benefit  of  post-test  analysis.  It  is 
very  important  that  the  range  meteorologist  interact  with  the  project  directors 
early  in  the  mission  planning  stage  to  define  meaningful  and  achievable  measure¬ 
ment  goals  and  to  identify  mission  meteorological  limitations.  Regulatory  agencies 
are  also  showing  increasing  interest  in  measurement  relevance,  particularly  m  the 
air  quality  area.  Hoffhagle  et  al.  (1981),  for  ex^ple,  de^e  relevant  instrumenta¬ 
tion  and  information  requirements  to  characterize  diffusion  fi’om  point  sources. 

14  2  Site  Selection  and  Instrument  Exposure.  Site  selection  is  an  important 
aspect  of  relevance  for  meteorological  measurements.  To  provide  representative 
data,  instruments  must  be  mounted  in  an  environment  where  they  can  perform 
according  to  their  design  specifications.  Functional  testing  can  help  determme 
which  instruments  can,  in  a  given  setting,  provide  the  most  representative  date. 
The  site  chosen  for  measurement  must  adequately  expose  the  instn^ent  to  the 
phenomenon  to  be  measured.  The  site  must  also  be  one  where  tee  instrument  is 
not  likely  to  suffer  damage  or  tampering.  Site  selection  in  flat,  isolated,  open 
terrain  is  not  difficult,  but  difficulties  increase  in  tee  presence  of  complex  teirain 
and  manmade  obstacles.  If  sensor  exposure  is  compromised  by  other  operational 
constraints,  a  description  of  the  problem  and  its  anticipated  effects  should 
accompany  date  obtained  fi'om  tee  site. 

Instrument  exposure  can  be  defined  in  terms  of  fetch,  which  is  the  distance 
upwind  fi-om  a  measurement  site  to  tee  nearest  discontinuity  or  obstacle  that  is 
likely  to  influence  instrument  readings.  Wind  measurements  are  tee  most 
susceptible  to  exposure  problems.  For  wind  readings,  a  fetch  to  obstacle  height 
ratio  of  at  least  15:1  is  needed  to  obtain  undisturbed  wind  profiles  (Wieringa, 
1976).  The  measurement  site  should  also  be  at  least  150  meters  fi*om  any 
discontinuity  in  site  roughness.  To  obtain  representative  measurements  above  a 
surface,  an  anemometer  should  be  mounted  at  a  height  at  least  five  times  tee 
surface  roughness  length  above  the  surface.  The  long-term  correlation  between 
wind  speed  and  wind  direction  should  be  at  its  maximum  along  tee  prevailing 
wind  direction  for  a  properly  exposed  anemometer  stationed  in  flat,  open  terrain. 
A  persistent  correlation  between  wind  speed  and  direction  that  is  imexplained  by 
the  prevailing  winds  can  provide  a  quantitative  measure  of  tee  terrain  or  obstacle 
influence  on  wind  measurements.  Quantitative  measures  of  persistence  and 
correlation  are  described  in  chapter  3,  paragraph  3.2.  A  detailed  discussion  of 
roughness,  fetch,  and  rules  for  matching  a  measurement  array  to  the  available 
terrain  is  presented  by  Wieringa  (1993).  Equipment  siting  guidance  can  be 
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obtained  from  the  manufacturer,  from  a  command  techmcal  assistance  office,  or 
from  a  standards-writing  organization  such  as  the  American  Society  for  Testing 
and  Materials  (ASTM). 

2  4  3  gfl-mnling  and  Averaging.  To  produce  relevant  information,  measurements 
must  be  made  using  instruments  with  suitable  response  characteristics  and 
appropriate  sampling  rates.  Instrument  response  is  usually  defined  in  terms  of  a 
time  or  distance  constant.  This  response  characteristic  should  exceed,  by  a  factor 
of  10  or  greater,  the  scale  or  duration  of  the  least  significant  event  to  be  ^ 
measured.  The  sampling  rate  should  be  comparable  with  the  instrument’s 
response.  For  example,  there  is  little  to  be  gained  by  sampling  once  per  second  or 
trying  to  resolve  features  of  less  than  2-minutes  duration  using  a  thermometer 
with  a  12-second  response  time.  Choosing  an  instrument  with  response  charac¬ 
teristics  and  a  sampling  rate  consistent  with  the  anticipated  measurement 
requirements  maximizes  the  available  information  content  while  minimizing 
measurement  redundancy. 

Averaging  time  should  be  a  function  of  application  and  atmospheric 
conditions.  Data  collected  within  a  specified  averaging  time  are  used  to  produce 
statistics;  these  statistics  constitute  the  information  needed  to  evaluate  the  data 
and  model  the  phenomenon  of  interest.  Averagmg  should  be  avoided  over  periods 
when  major  trends  or  step  changes  are  present  in  the  data. 

The  choice  of  an  averaging  period  is  also  dependent  on  sampling  rate  and 
the  statistics  of  interest.  First-moment  statistics  (means)  require  fewer  data 
points  to  produce  reliable  estimates  than  do  second-moment  statistics  (variances). 
For  example,  given  an  anemometer  and  vane  sampling  at  a  rate  of  1  Hz,  1-minute 
averaging  may  be  sufficient  to  define  a  mean  wind,  but  15  minutes  of  data  may  be 
required  to  define  wind-direction  variance  with  an  acceptable  level  of  statistical 
significance.  Lenschow  et  al.  (1994)  provide  guidance  on  choosing  suitable 
averaging  times  and  present  a  methodology  for  defining  random  versus  systematic 
errors  in  turbulence  statistics. 

While  sufficiently  long  averaging  periods  are  needed,  excessively  long 
averaging  times  and  large  sample  sizes  should  be  avoided.  Wind  direction 
averaging  in  excess  of  1  hour  may  include  undesired  long-term  trends.  The  effects 
of  trends  can  be  minimized  through  the  use  of  detrending  (high-pass  filter) 
algorithms.  Detrending  is  particularly  important  for  spectrum  analysis  (see 
chapter  4,  paragraph  4.5).  Conversely,  unnecessary  high  frequency  data  or 
random  noise  can  be  suppressed  through  block  averaging  or  low-pass  filters. 

Block  averaging  produces  independent  samples,  while  running-mean  and  median 
filters  provide  smoother  point-to-point  transitions. 
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14  4  Information  Dissemination.  Computers  are  invaluable  mds  in  ^e  date 
collection  analysis,  display,  and  archival  processes,  but  with  this  ad^tional  speed 

Smes  the  e=cpe^tion  that  date  wiU  n^e  avaalable 
ouicklv.  Date  should  not  be  disseminated  without  quahty  control.  Raind 
dissemination  of  dubious  or  irrelevant  date  is  not  in  the  best  interest  of  range 
“erf  A  statement  of  quality  and  an  explanation  of  QC  P-^edures  ^ 
should  accompany  all  disseminated  information.  A  review  by  a  meteorologist  with 
knowledge  of  mission  objectives  and  information  requirements  and  a  firm  under- 
standing  of  meteorological  processes  remains  the  key  to  a  successful  measurement 

program. 

Date  that  are  meant  to  be  shared  with  other  users  should  be  presented,  if 
possible,  in  ASCH  format.  A  multitude  of  media  are  available  for  ardm^g  date, 
but  a  9-track  tape  remains  the  only  medium  that  is  common  to  ^ 
for  the  transfer  of  large  quantities  (tens  of  megabytes  or  more)  of  date,  although 
compact  discs  (CD-ROM)  are  emerging  as  a  common  mass-storage  metei^. 
Apphcations  of  CD-ROM  in  the  atmospheric  sciences  are  discussed  by  Mass 
(1993)  For  smaller  quantities  of  data,  floppy  disks  are  inexpensive  and 
convenient  and  can  be  read  at  most  ranges.  Data  compression  enh^ces  the 
storage  capacity  of  floppy  disks.  The  3^-inch  disks  are  more  robust  and  generally 

survive  transportation  better  than  5V4-inch  disks. 

Date  formats  must  include  a  certain  minimum  amoimt  of  docu^ntary 
information  to  be  of  value  to  anyone  other  than  the  immediate  user.  This 
rniniTmiTTi  information  includes  station  identification,  a  date/time  stamp,  the 
sampling  rate  and  averaging  interval,  sensor  height  above  ^oimd,  sensor  location, 
and  measure  units.  Data  sets  intended  for  future  use  should  also  mdude 
documentation  on  instrument  type,  axis  orientation  (as  applicable),  ^te 
corrections  or  adjustments,  sampling  interval,  project  name,  md  date  quahly 
indicators.  Additional  test  information  should  mclude  a  site  descnpUon  conmsfing 
of  roughness  elements,  obstacles,  major  terrain  features,  and  the  probable  effects 
of  these  features  on  the  date  such  as  wake  effects  for  certam  wmd  directions. 

Date  presentation  should  be  dear  and  in  standard  meter-l^ogram-second  (mks) 
units  unless  the  use  of  alternative  units  is  coordinated  with  customers.  S^d^d 
mks  units  are  presented  in  Practice  for  Use  of  the  International  System  of  Unite 
fST^  ftbo  Modernized  Metric  Svsteml  (ASTM  E380).  Time  is  preferably  presented 
in  Universal  Coordinated  Time  (UTC).  Local  standard  time  (LST)  is  acceptable  for 
date  sets  designed  primarily  for  local  use.  Time  stamps  should  mdicate  whether 
UTC  or  LST  is  used.  For  data  averaged  over  time,  the  time  sl^p  is  usually 
taken  at  the  beginning  of  the  averaging  period.  If  an  alternative  convention  is 
used,  that  convention  should  be  noted  in  the  date  set. 
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Hard-coDV  data  sxxmmaries  are  more  meaningful  when  presented  ^  ® 
format  designed  for  user  interpretation.  For  example,  wind  speeds  and  dire<*ons 
^a  Station  are  most  useful  when  presented  in  a«a,^nt  eolu^  and  m 
7Zly  dLfined,  common  units.  The  same  is  usually  true  for  temperatures  and 
dew  points.  Exceptions  are  made  for  some  appUcations.  On  a  mulUevel^er. 
^d^eeds  and  temperatures  may  be  best  presented  m  an  adjacent 
format  when  the  gradients  of  these  variables  me  of 

formats  are  most  easily  read  when  there  is  sufficient  separataon  (at  least  two 
snaces)  between  columns.  The  number  of  significant  digits  should  not  exceed  ffie 
resolution  of  measurements.  For  example,  wind  speed  and  t^peratare  ^  Iffidy 
to  be  significant  to  a  tenth  of  a  unit,  while  wind  direcbon  IJ^ly  be 
to  the  nearest  whole  degree.  Obtaining  customer  input  on  desired  presentataon 
formats  before  measurements  are  taken  can  save  time  and  improve  customer 


satisfaction. 
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CHAPTER  2 


ERROR  SOURCES 


2.1  Missing  Data 

Because  an  "open  line"  or  "short"  can  sometimes  be  confused  with  a  valid 
zero  or  full-scale  voltage  reading  when  using  analog  instruments,  a  method  to 
uniquely  identify  the  condition  when  no  data  are  present  is  a  concern.  Zero  or 
full-scale  voltages  are  easier  to  identify  before  conversion  to  their  equivalent  m 
engineering  units.  Likewise,  a  single  bad-data  point  is  usually  easier  to  identify 
before  it  is  averaged  with  other  data  points.  It  is,  therefore,  desirable  to  perform 
checks  for  open  hne  or  short  conditions  early  in  the  QC  process  before  averagmg 
or  conversion  to  engineering  units.  For  anemometers,  0  to  5  volts  mght  convert  to 
a  speed  range  of  0.2  to  50  meters  per  second  (m/s).  Mtematively,  if  the  voltage 
range  is  -5  to  +5  volts,  zero  volts  may  produce  a  reading  of  25  m/s.  For  a 
thermometer,  0  volt  is  often  set  for  the  lowest  temperature  of  the  operatmg  range, 
producing  a  reading  of  -49.9  °C,  while  full  scale  is  equivalent  to  50  °C.  Wind 
direction  can  be  particularly  confusing  because  either  0  or  full  scale  (4.9  volts) 
might  indicate  a  north  wind  reading  from  a  sensor  with  a  single  potentiometer. 
Wind  direction  monitors  using  dual  potentiometers  often  set  the  first  poten¬ 
tiometer  to  range  between  0  and  270°  and  the  second  potentiometer  to  range 
between  270  and  540°.  Identifying  the  "no  data"  condition  is  a  difficult  problem, 
particularly  for  an  automated  data-collection  system. 

Given  the  difficulty  of  identifying  missing  data  versus  a  valid  null  voltage 
reading,  analysis  of  multiple  sensor  records  can  help  identify  a  missing  data 
condition.  For  example,  a  typical  mesometeorological  network  (mesomet)  station 
might  report  wind  speed,  wind  direction,  and  temperature.  A  0  volt  reading  m^ 
be  valid  data  for  one  sensor  and  may  even  persist  for  a  few  hours  (calm  winds,  for 
example).  However,  it  is  unlikely  that  the  atmosphere  would  produce  conditions 
that  register  as  null  voltages  simultaneously  for  two  or  more  sensors  on  a  given 
station;  that  is,  a  persistent  north  wind  of  0.2  m/s  or  less  and  a  temperature  of 
-49.9  °C  is  unlikely  even  if  the  station  is  located  in  the  Yukon.  For  most  locations, 
the  simultaneous  acquisition  of  nuU-voltage  data  points  on  two  or  more  sensors  at 
a  single  mesomet  site  provides  strong  evidence  of  a  no-data  "open  line"  condition, 
particularly  if  these  readings  persist  over  several  averaging  intervals.  Analogous 
criteria  can  be  used  to  define  a  "short"  condition.  Monitoring  an  on-site  battery 
reference  voltage  can  be  helpful  in  identifying  the  existence  of  a  no-data  condition 
at  the  cost  of  logging  an  additional  channel  of  data. 
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Line  noise  may  also  be  present  to  produce  some  dither  (irreducible  r^dom 
noise)  or  bias  voltage  even  when  nothing  is  coming  from  a  sensor.  Thus,  hne  noise 
may  produce  the  appearance  of  sensor  activity  even  when  an  open  hne  or  short 
condition  exists.  Consequently,  it  is  important  to  periodically  perform  zero-voltege 
checks  through  the  system.  If  a  dither  of,  say,  ±0.01  volt  or  a  bias  of  +0.02  volt  is 
found,  it  should  be  considered  in  the  QC  criteria  for  null  data. 

Automated  procedures  should  be  used  to  flag  suspected  null  data.  If  this  ^ 
condition  is  confirmed  by  a  meteorologist  during  the  QC  process,  these  data  points 
should  be  assigned  a  unique  character  or  numerical  value  that  is  easily  disfrn- 
guished  within  the  archived  data  set.  The  preferred  identifier  for  missing  data  is 
a  string  of  9999.  The  string  should  fill  all  the  digits  available  for  any  data  charac¬ 
ter  (that  is,  -49.9  becomes  999.9).  Alternatively,  a  special  control  character  ^  be 
used  to  flag  suspect  data  points.  The  #  sign  can  be  used  for  this  purpose  (-49.9 
becomes  -49.9#).  The  advantage  of  using  a  control  symbol  of  this  sort  is  that  the 
analyst  can  observe  each  reported  value  and  add  or  remove  control  symbols  as 
deemed  necessary;  data-processing  routines  can  be  programmed  to  ignore 
numerical  values  followed  by  a  #  sign. 

An  open  line  or  short  condition  is  less  likely  to  be  confused  with  real  data  for 
digital  instruments  because  the  transmitted  character  is  a  string  of  0  or  1.  Digital 
transmission  is  also  done  in  a  specific  message  format  that  precludes  infusion 
with  an  open  line  or  short  condition.  Received  data  packets  can  be  subjected  to 
format  (check  sum)  and  parity  checks  that  unambiguously  identify  gross-error 
conditions.  In  this  respect,  some  basic  QC  is  already  built  into  digital 
instrumentation.  On  the  other  hand,  the  digitally  transmitted  signal  may  mdude 
pre-averaged  data  in  which  noise  or  short  conditions  are  already  averaged  mto  the 
transmitted  signal,  creating  a  very  dififlcult  QC  problem. 

2.2  T^nisy  nomnuinications  Lines  or  Acquisition  System 

Analog  output-voltage  signals  are  susceptible  to  contamination  during 
transmission  from  the  sensor  to  the  point  where  data  averaging  and  storage  are 
performed.  Analog  data  typically  originate  as  voltage  or  current  signals  from  the 
sensor  head  arid  pass  along  a  wire  to  a  translator  that  amplifies  and  conditions 
the  signal,  which  is  then  passed  to  a  data-logging  device. 

Communications  fines,  particularly  metal  cables,  present  many  opportunities 
for  signal  corruption.  Cable  resistivity  causes  a  voltage  loss,  imposing  a  bias  on 
the  signal.  All  electrical  systems  are  referenced  to  ground  and  are  susceptible  to 
grounding  changes  because  of  variations  in  soil  moisture  content.  Moisture 
penetration  into  the  cable  sheath  may  also  impose  variable  voltage  offsets.  Cables 
also  act  as  antennas,  receiving  inductance  from  other  nearby  cables,  p^icularly 
power  cables,  and  are  susceptible  to  changes  in  the  atmospheric  electric  field. 
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Some  of  these  conditions  can  be  identified  and  corrected  through  end-to-end  si^al 
checks,  but  others  are  variable  and  must  be  minimized  through  careful  groxmding 

and  cabling  practices. 

Good  grounding  and  cabling  practices  are  necessary  for  the  acquisition  of 
analog  data  fi-ee  of  noise  and  bias.  Signal  cables  should  never  He  paraUel  to  power 
cables  or  in  long  runs  next  to  other  signal  cables,  particularly  when  they  are 
poorly  shielded.  When  it  is  necessary  to  have  cables  near  each  other,  they  should 
be  shielded.  Cabling  should  be  laid  out  to  cross  at  right  angles  to  minimize 
inductance.  Also,  common  and  stable  grounding  should  be  established  for  all 
equipment  on  a  tower.  If  possible,  this  groxmding  should  be  common  with  the 
data  acquisition  system.  If  common  grounding  is  not  achieved,  ground  loop 
voltages  may  pass  through  the  system.  Grounding  and  ground  loop  conditions 
often  vary  with  changes  in  soil  moisture. 

Translators,  signal  conditioners,  and  ampHfiers  contribute  a  certain  amount 
of  background  noise  and  bias  to  an  analog  signal.  The  signal  at  this  point  is  also 
susceptible  to  field  effects  because  of  radio  transmitters  or  other  ne^by  sources. 
Radio  transmitters,  when  keyed,  can  produce  data  spikes.  Electronic  interference 
checks  shotdd  be  done  by  systematically  turning  on  or  ke^g  equipment  while 
watching  for  signal  output  level  changes  in  the  data  monitoring  system. 

Signals  transmitted  as  frequency  coimts  or  in  digital  form  are  less  subject  to 
noise  and  interference  than  are  voltages.  Wind  speed  signals,  often  transmitted  as 
pxilses,  are  less  susceptible  to  errors  caused  by  groxmding  and  cabling  faxilts  than 
wind  direction  voltages.  Digital  signal  qxiality  is  also  insensitive  to  small  changes 
in  voltage  level  and  is  less  susceptible  to  these  errors. 

2.3  Environmental  Hazards 

In  addition  to  electromagnetic  hazards,  the  biological  and  meteorological 
environment  creates  special  hazards  for  meteorological  eqxiipment,  particxilarly  for 
eqxiipment  left  in  the  field  for  extended  periods.  Insects  bxiild  nests  or  webs  in 
exposed  hoxisings,  air  intake  ducts,  sxmshades,  and  electronics  shelters 
openings.  The  first  indication  of  a  problem  may  come  as  an  intermittent  open  or 
"short"  in  the  signal,  but  the  instrument  can  experience  slow  degradation  of  signal 
quality  well  before  a  faxilt  is  detected.  Rodents  are  known  to  chew  cables  and 
bxirrow  xmder  or  into  standing  shelters.  Birds  perch  on  xvind  vanes  bending  the 
arms  and  attack  rotating  instrximents  such  as  cup  anemometers. 

Sites  with  known  biological  hazards  will  reqxiire  special  attention  such  as 
seaHng  electronics  housings  and  designing  eqxiipment  moxmts  to  make  animal 
access  difficult  or  xmattractive.  Installing  a  bird  perch  above  xvind  eqxiipment  can 
reduce  the  incidence  of  bent  wind  vanes.  Cup  anemometers  that  strongly  reflect 
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sunlight  seem  to  provoke  attack  by  raptors,  and  many  insects  are  attracted  to  the 
color  yellow.  Choosing  colors,  coatings,  and  configurations  that  do  not  provoke  or 
attract  birds  or  insects  can  help  minimize  damage.  Cabling  at  a  permanent  site 
can  be  buried,  although  the  caution  against  long  runs  of  closely  Imd  p^allel  cable 
remains.  Cabling  lying  on  the  ground  is  susceptible  to  gnawing  by  rodents. 
Mounting  cables  on  poles  above  the  ground  may  provide  adequate  cable  separa- 
tion,  but  cable  coverings  should  be  designed  to  withstand  damage  caused  by  wmd 
loading  and  degradation  by  ultraviolet  radiation. 

Meteorological  effects  include  sensor  degradation  caused  by  intrusion  or 
accumulation  of  moisture,  salt,  dust,  mud,  or  ice  and  d^age  fi-om  overheatmg  or 
severe  weather.  Rotating  instruments  such  as  mechamcal  anemometers  and  wmd 
vanes  are  particularly  susceptible  to  these  effects.  Salt  accumulation  is  prevalent 
in  the  marine  environment  where  sea  spray  leaves  a  deposit  on  instr^ents.  In 
dry  environments,  fine  dust  can  filter  into  electronics  wmponente  and  bearing 
assembbes,  and  long-term  exposure  to  intense  ultraviolet  radiation  causes  ptostic 
or  composite  components  to  become  brittle  and  crack.  Mud  splatter  during  heavy 
precipitation  events  can  coat  instruments  and  block  aspirators  or  cooling  vents. 

Ice  in  the  form  of  rime  or  snow  can  accumulate  on  sensors,  and  ice  m  the  form  o 
hail  and  strong  winds  can  physically  damage  instrumentation. 


Humidity  is  a  major  environmental  hazard  for  instrumentation.  Condensa¬ 
tion  in  high  humidity  conditions  can  cause  electronics  faults  or  distortion  on 
optical  surfaces.  Mold  and  fimgi  can  also  grow  on  damp  surfaras,  degrading  si^al 
strength  or  reducing  sensitivity.  Growth  often  begins  on  shields  or  shelters  and 
may  go  imdetected  until  noticeable  instrument  performance  degra^tmn  occurs. 
Very  low  humidity  is  accompanied  by  increased  static  electricity  which  ^ 
damage  electronics  or  cause  dust  to  accumulate  on  sensors  or  optical  surfaces. 

Environmental  hazards  impact  the  quality  of  meteorological  data  in  a  variety 
of  ways,  depending  on  the  sensor.  Cup  or  propeller  anemometers  may  cease  to 
turn  after  receiving  deposited  material,  or  may  break  fi:ee  of  the  material  only 
with  the  passage  of  a  sudden  gust  of  wind.  Accumulations  of  foreign  ^tter  can 
also  alter  an  instrument’s  aerodynamic  characteristics,  thereby  degrading  its 
threshold  and  response  distance.  On  nonmechanical  instruments,  these  hazards 
can  appear  as  signal  loss,  bias,  or  noise  because  of  ventilation  duct  blockage  or 
obscuration  of  the  path  between  an  optical  transmitter/receiver  pair. 


Adverse  meteorological  effects  on  data  quality  can  be  mii^zed  through 
careful  selection  of  equipment  designed  for  use  in  adverse  environmente.  Ice 
accumulation  can  be  minimized  by  using  heater  elements  and  by  covering  sensor 
surfaces  with  a  hydrophobic  coating.  Frequent  cleaning  for  sea  spray  or  dust  and 
special  site  visits  after  storms  may  be  necessary  to  keep  mstruments  m  good 
operating  condition. 
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2.4  Aliimmcnt  Errors 

Alignment  accuracy  is  particularly  important  for  wind-measurement  equip¬ 
ment  Anemometers  not  vertically  aligned  will  exhibit  unusual  acceleration 
characteristics  (MacCready,  1966).  A  vrind-vane  aUpiment  error  appeare  M  a  bias 
in  the  wind  direction  measurements.  Mechamcal-wmd  equipment  sho^d  1« 
aligned  to  within  ±2°  of  vertical  to  minimize  measurement  hip  from 
winds.  Instruments  set  up  to  measure  alonpvind  (u),  crossed  (v),  ^d  vertical 
wind  (w)  components  also  require  careful  alignment  (to  withm  ±0.1  )  for  stress 
measurements  (Kaimal  and  Haugen,  1969).  Cross-component  contamination 
caused  by  misalignment  is  difficult  to  remove.  Vertic^  velocity,  usuaUy  a  smaU 
fraction  of  the  horizontal  wind,  is  particularly  susceptible  to  cross-component 

contamination. 

Alignment  with  respect  to  either  magnetic  norffi  or  true  north  is  acreptable, 
but  this  alignment  choice  must  be  clearly  stated.  Wind  direction  is  usual  y 
reported  with  respect  to  true  north,  so  instruments  aligned  with  respect  to 
magnetic  north  wiU  need  to  have  a  software  correction  apphed  to  the  ^gular 
information  unless  the  magnetic  declination  is  zero.  Magnetic  dechnati^ 
information  presented  as  a  diagram  of  the  relationship  between  gnd  north,  true 
north,  and  magnetic  north,  can  be  found  on  U.S.  Geological  Su^ey  or  Defense 
Mapping  Agency  maps.  The  decUnation  angle  varies  slowly  vnth  time  bemuse  of 
changes  in  the  magnetic  pole,  so  the  use  of  fairly  recent  maps  (vnthin  the  last  20 
years)  and  the  yearly  correction  (also  annotated  on  these  maps)  is  recomnaended. 
If  the  magnetic  north  pointer  is  east  (west)  of  true  north,  add  (subtract)  the 
magnetic  declination  to  (from)  the  magnetic  north  wind  reading  to  get  a  wind 
reading  with  respect  to  true  north.  AHgnment  aci^acy  should  be  mthm  ±2  . 
Lockhart  (1989A)  provides  a  methodology  for  obtaining  a  precise  ahgnment  usmg 

the  true  solar-noon  method. 


2.6  r.nlihration  Inadequacies 

Static  and  dynamic  calibrations  are  used  on  instruments  with  analog  output. 
Static  calibration  usually  involves  mounting  an  instrument  in  a  test  fixture  and 
checking  voltage  output  at  a  series  of  predetermined  positions.  Static  cahbrations 
check  the  validity  of  algorithms  used  to  convert  voltages  to  engmeenng  umts  orer 
the  operating  span  but  provide  no  information  on  the  dynamic  perform^ce  of  the 
instrument.  Static  cahbration  is  adequate  for  instruments  measui^g  thermo- 
dynamic  quantities  such  as  temperature,  pressi^e,  or  humidity  but  is  msi^aent 
for  the  calibration  of  rotating  wind  speed  and  direction  eqmpment.  ^t^ng-wmd 
equipment  require  both  mechanical  and  electrical  checks.  When  cahbrations  are 
done,  electrical  circuitry  designed  to  dampen  signals  (usually  resistor  rapaator 
(RC)  filters)  should  be  bypassed  so  that  the  true  electromechamcal  mstrument 

response  is  measured. 
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Dynamic  caHbration  should  be  performed  in  a  wind  tunnel  where  anemom¬ 
eter  and  wind-vane  performance  can  be  defined  with  respert  to  steady-state  flow 
conditions.  Dynamic  measures  include  the  threshold,  which  is  the  lowest  p)eed  at 
which  a  rotating  anemometer  starts  and  continues  to  hm,  or  a  vane  wtoch  s^ 
to  turn  towards  the  true  wind  direction  fi’om  an  initial  displacement  of  10  .  Ihe 
threshold  should  be  checked  periodically  by  the  user,  because  it  is  a  sensitive 
indicator  of  bearing  degradation.  Other  dynamic  measures  fo™d  v^es  mdude 
delay  distance,  damping  ratio,  and  overshoot  as  defined  by  AS™-D5096  (1990) 
and  Lockhart  (1989A).  Performance  and  calibration  specifications  should  be 
provided  by  the  manufacturer. 


Although  the  dynamic  performance  of  wind  equipment  is  defined  in  a  low- 
turbulence  wind  tunnel,  its  actual  performance  in  a  turbulent-wmd  field  is 
another  matter.  Mechanical  wind  equipment  is  cahbrated  after  haymg  the 
opportunity  to  come  into  equilibrium  with  the  tunnel’s  velocity  field  If  subjected 
to  atmospheric  turbulence  at  scales  comparable  to  or  greater  than  ihe  distance 
constant  (L),  the  sensor  cannot  approach  equilibrium  and  will  provide  an 
attenuated  wind  reading.  A  rotating  anemometer  exposed  to  a  step  change  m 
wind  speed  (Au)  over  a  time  period  (t)  will  report  an  attenuated  fi-action  (A)  ot  the 

Au  expressed  by 


A  =  [1  +  ilnULutfr^ 


(2-1) 


Therefore,  wind  equipment  ivith  a  large-distance  constant,  even  if  well  calibrated, 
cannot  provide  representative  wind  measurements  near  the  surface  or  m  locahons 
where  gusty  conditions  are  prevalent.  Near-surface  wind  measurements  require 
fast-response  sensors,  while  heavier  slow-response  sensors  are  more  appropriate 
for  wind  measurements  at  greater  heights  where  stronger  but  less  turbulent  wmd 

conditions  occur. 

In  addition  to  providing  current  calibrations  for  equipment,  cahbration 
fadhties  should  report  the  calibration  errors  found  when  an  instrument  entered 
the  facility  These  errors  can  then  be  used  to  correct  data  collected  by  the 
equipment  prior  to  recalibration.  Also,  the  record  of  cahbration  errors  can  be  used 
to  define  a  cahbration  schedule.  Instruments  that  consistently  extobit  no 
cahbration  drift  can  be  put  on  a  longer  cahbration  interval,  while  instruments  that 
consistently  exhibit  significant  calibration  drift  should  be  returned  to  the  cahbra¬ 
tion  fadlity  sooner  (MIL-STD  45622A,  1980). 
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Time  synchronization  is  necessary  for  all  meteorological  measurements,  and 
time  should  be  part  of  any  measurement  record.  Time-accuracy  requirements  vary 
with  the  application.  Accuracy  to  within  a  minute  is  adequate  for  synoptic 
appUcations,  but  accuracy  to  within  a  tenth  of  a  second  is  needed  if  the  recorded 
data  are  to  be  used  for  flux  computations.  Computer-system  time  is  often  derived 
from  the  60-Hz  signal  passing  through  the  power  supply.  A  more  adequate  time 
base  is  obtained  by  reference  to  IRIG  timing  or  to  a  GPS-satellite  time  base.  The 
reference  time  system  being  used  should  be  identified.  Universal  Coordinated 
Time  (UTC)  is  the  preferred  dock  time,  espedally  for  data  that  are  to  be  iised  off 
range.  If  UTC  is  not  chosen  for  on-range  applications,  local-standard  time  is 
recommended  for  use  throughout  the  year. 
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CHAPTER  3 


REAL-TIME  QUALITY  CONTROL 


3.1  OperatioTinl  Quality  Control  Checks 

Operational  qiiality  control  checks  to  monitor  data  qu^ty  include  those 
actions  that  be  taken  during  or  shortly  after  data  collection.  These  actions 
include  using  one’s  knowledge  of  the  measurement  system  and  the  measured 
geophysical  phenomena  during  inspection  of  the  data  and  summa^  statistics,  and 
using  real-time  predictors  to  flag  possible  fault  conditions.  Steps  in  the  QC 
process  are  outlined  briefly  below. 

3.1.1  Become  familiar  with  instrument  and  data-system  performance  character¬ 
istics,  error  modes,  and  tolerances.  Much  data  are  lost  because  the  operator  does 
not  imderstand  equipment  capabilities  and  limitations.  Each  instrument  has 
characteristic  error  modes  that  a  well-trained  analyst  should  be  able  to  identify 
during  a  review  of  the  data. 

3.1.2  Perform  out-of-range  and  gross  departure  checks  as  dose  to  the  measme- 
ment  source  as  possible.  Errors  masked  by  averaging  or  conversion  into  derived 
variables  are  more  diflBcult  to  detect  and  cause  more  information  loss  than  errors 
detected  early  in  the  measurement  process. 

3.1.3  Use  known  physical  relationships  to  check  measurement  validity.  Flag 
gross  departures  ftom  the  hydrostatic  condition,  the  adiabatic  lapse  rate,  or  large 
step  changes  in  signal  level.  These  conditions  represent  either  an  emerging  fault 
or  the  onset  of  an  interesting  meteorological  condition.  In  either  case,  these 
conditions  warrant  further  attention. 

3.1.4  Use  intersensor  and  intervariable  comparisons.  Changes  in  one  measured 
field  are  often  accompanied  by  changes  in  another.  For  example,  a  change  in 
pressure  should  be  accompanied  by  a  change  in  the  wind  field.  A  major  change  in 
one  variable  without  an  accompanying  change  in  related  variables  may  be  due  to 
the  emergence  of  an  error  or  fault  condition. 

3.1.5  Plot  time-series  of  the  data  for  visual  inspection.  Plotted  data  c^  be 
quickly  scanned  for  trends,  spikes,  or  flat  response  conditions.  If  possible,  plots 
should  encompass  one  or  more  meaningful  cycles  such  as  the  diurnal  heating 
cycle.  Profile  plots  can  be  quickly  scanned  for  faults.  For  example,  figure  3.1 
contains  a  profile  of  four  unedited  v-component  (crosswind  axis)  sonic  anemometer 
data  sets  collected  simultaneously  at  a  rate  of  5  Hz  over  a  period  of  1  hour  (18,000 
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data  points  in  each  set).  The  top  panel  is  the  v-component  meas^ed  at  the 
2-ineter  level  on  a  meteorological  tower,  while  the  lower  panels  show  ^e  v- 
components  measured  respectively  at  4-,  6-,  and  8-meter  on  the  same  t^er^ 
quick  visual  inspection  reveals  clusters  of  noise  spikes  at  the  4-  and  8-  meter 
levels,  and  two  isolated  spikes  on  the  6-meter  level. 

3.1.6  Calculate  summary  statistics  over  meaningful  averaging  periods.  Choice  of 
a  meaningful  averaging  period  depends  on  the  time  scale  of  phenomenon  o 
interest  but  should  include  at  least  250  data  points  for  escalation  of  a  stan^d 
deviation.  The  basic  summary  statistics  (the  mean,  median,  ma^um,  mimm^, 
and  standard  deviation)  calculated  for  this  averaging  period  can  be  good  indicators 
of  the  presence  (or  absence)  of  noise. 


3  17  Use  the  last  measured  data  point  and  a  running  mean  of  10  to  100  data 
points  as  a  real-time  predictor  for  each  new  data  point.  Signific^t  deviations 
from  the  predicted  value  can  indicate  either  the  emergence  of  a  fault  condition  or 
a  change  in  meteorological  conditions. 


3.2  Meteorologies*!  Variables  Basic  Statistics 

Meteorological  variables  consist  of  scalar  quantities  such  as  temperature  and 
humidity  and  vector  quantities  representing  the  wind.  Basic  statistics  for  sc^ar 
quantities  can  be  computed  using  methods  presented  in  statistics  texts.  Wmd 
data  can  be  presented  as  components  (eastings  and  northings)  or  as  speed  and 
direction.  The  treatment  of  wind  requires  special  care  because  of  the  discontinuity 
in  the  wind  direction  circular  scale  between  360°  and  001°  and  because  of  possible 
ambiguity  in  the  quadrant  representation  of  the  arctangent  ^ction.  Also,  the 
choice  of  algorithm  (for  example,  scalar  or  unit  vector)  will  bias  the  outcome. 
Algorithm  selection  should  be  based  on  how  the  wind  data  are  to  be^ed.  A 
scalar-averaged  wind  speed,  unit  vector-averaged  wind  direction,  and  ^cul^ 
standard  deviation  are  appropriate  for  most  appHcations.  The  rem^der  of 
paragraph  3.2  contains  computationally  efficient  single-pass  algorithms  for 
computation  of  scalar  or  vector  quantities.  Table  3-1  shows  sonic  anemometer- 
derived  summaiy  statistics  obtained  by  using  many  of  these  algorithms. 


3.2.1  Scalar  Mean.  For  n  samples  of  a  scalar  quantity  a,  the  mean  (o)  is  the 
sum  of  the  samples  divided  by  n. 


(3-1) 


3-3 


If  the  quantity  to  be  computed  is  a  scalar  mean  horizontal  wind  speed  (w) 
obtained  from  along-axis  (u-component)  and  cross-ass  (v-component)  measure- 
ments,  the  mean  is  the  sum  of  the  squares  of  each  individual  velocity  component. 
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3  2  2  Vprtnr  Mean  Speed.  The  horizontal  vector  mean  wind  speed  (wv)  is  the_ 
square  root  of  the  sum  of  the  squares  of  the  mean  alongaxis  u  and  crossaxis  M 
horizontal  velocity  components. 


n 


(3-3) 


NOTE:  The  vector  mean  wind  speed  is  applicable  to  calculation  of  atmospheric 
transport  for  diffusing  clouds.  The  scalar  mean  wind  speed  is  the  variable  of 
choice  for  most  other  applications.  A  ratio  of  the  vector  to  scalar  wmd  speeds  is  a 
measure  of  wind  direction  persistence  (Panofsky  and  Brier,  1965). 

3.2.3  Unit  Vector  Moan  Direction.  For  a  set  of  wind  direction  angle  (Oj)  observa¬ 
tions,  the  unit  vector  mean  wind  direction  (0)  is  obtained  from  the  arctangent  of 
the  averaged  sines  (E)  and  cosines  (N)  of  the  wind  direction  angles 


£  =  —  V  sine; 

n  U 


—  1  ^ 

N  =  —  y)  cosS; 
n  i=i 


(3-4) 


(3-5) 


0'  =  arctan  (N/JE) 


(3-6) 
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To  ensure  that  the  resultant  falls  within  the  correct  wind  direction  quadrant, 


e  = 


I  90  -  e',  E<o 
270  -  e',  E>0 


(3-7) 


NOTE-  The  negative  signs  in  equations  (3-4)  and  (3-5)  are  used  to  achieve  the 
standard  meteorological  convention  with  wind  directions  from  the  west  and  south 
having  a  positive  sign. 

3  2.4  Sneed-Weighted  Mean  Direction.  The  wind  speed-weighted  mean  wind 
direction  includes  the  set  of  wind  speed  (Sj)  in  the  computation  of  sines  and 
cosines, 


S  5,  cosOi 
n  1=1 


smOi 

n  J.1 


and  processing  the  resulting  components  through  equations  (3-6)  and  (3-p. 
NOTE:  The  degree  to  which  the  unit  vector  and  speed-weighted  mean  directions 
agree  is  a  measure  of  the  lack  of  correlation  between  fluctuations  of  wind  speed 
and  direction  (Yamartino,  1984). 

3.2.5  Scalar  Standard  Deviation.  For  n  samples  of  a  scalar  quantity  a,  the 
standard  deviation  (Panofsky  and  Brier,  1965)  is 


o 


a 


(3-10) 


where  a  is  as  defined  in  equation  (3-1). 
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3  2  6  Persistence  Estimator  of  Standard  Deviation.  An  estimate  of  the  wind  ^gle 
standard  deviation  Oe  can  be  obtained  using  persistence  defeed  as  ^e  rafe  of  toe 
scalar  to  vector  mean  wind  speeds  as  defeed  by  equations  (3-2)  and  (3-3)  (Weber, 

1991). 


Og  =  105.75 


1  - 


'  wv\ 
^ws) 


1.5337 


(3-11) 


NOTE:  This  method  can  provide  a  Og  estimate  using  wind  component  data  when 
direct  wind  direction  measurements  are  unavailable. 


3.2.7  flirmlar  Standard  Deviation.  The  standard-deviation  estomtor  tfet 
produces  the  least  error  (±2“)  for  computation  of  the  circular  ™d  d^e^tion 
standard  deviation  (a,)  over  an  angular  range  of  0°  to  103.9“  (Yamartmo,  1984)  is 


Og  =  arcsin  (e)  [1  -  0.1547 


(3-12) 


where 

€  =  1  -  (£2  + 


(3-13) 


3.2.8  Covariance  and  Correlation.  A  covariance  (Cov)  is  toe  product  of  two  time- 
dependent  quantities.  For  quantities  A(t)  and  B(t)  with  means  of  A  and  B 


Cov  =  <[A{t)  -  A]  [B{t)  -  B}>, 


(3-14) 


where  <  >  indicates  a  user-selected  tone  average.  Covariance  normalized  by  toe 
variances  of  toe  time-dependent  quantities  gives  toe  coefficient  of 

correlation  (Cor)  between  these  quantities 


Cor 


<lA(t)  -  A\  mt)  -B]> 

[(oj)  (oj)^ 


(3-15) 
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Fiances,  standard  deviations,  tphbulence  intensities, 

AND  COVARIANCES, 


Date:  28  Feb  1994 
Trial  Nr.:  a69730 
Sonic  Loc.:  1.8  m  tri-axis 


Time:  0730  MST 
Location:  West  Tow 
Axis  Orient;  180"  (True) 


A.  MEAN  VALUES 


START 

TIME 

HR;MIN:SEC 

07:30:00: 

B.  VARIANCES 


SCALAR 

SPEED 

M/S 

0.4 


VECTOR 

SPEED 

M/S 

0.3 


WIND 

DIR 

DEG 

64.2 


ALONG 

AXIS 

M/S 

-0.125 


CROSS 

AXIS 

M/S 

0.259 


VERTICAL 

WIND 

M/S 

-0.02 


START 

UTP 

W 

WW 

TsTs’ 

REJECTED 

SPEED  OF 

TIME 

HR:MIN:SEC 

MW 

M^/S* 

M*/S* 

DEGK* 

DATA 

POINTS 

SOUND 

M/S 

07:30:00  0.112  0.055 

C.  STANDARD  DEVIATIONS  (Rotated) 

0.025 

0.044 

0 

336.6 

START 

SIGMA 

SIGMA 

SIGMA 

SIGMA 

SIGMA 

TIME 

U 

V 

W 

Ts 

THETA 

HR:MIN:SEC 

M/S 

M/S 

M/S 

DEG 

DEG 

07:30:00 

0.335 

0.235 

0.159 

0.211 

ND 

D.  TURBUT^ENCE  INTENSITIES  (RADIANS) 


START 

TIME 

HR:MIN:SEC 

07:30:00 


LONG 

Ix 

1.164 


LATERAL  VERTICAL 
ly  Iz 


E.  COVARIANCES  (Rotated) 


START 

TIME 

HR:MIN:SEC 

07:30:00 


inr 

MW 

-0.0005 


0.818 


UW 

MW 

-0.0171 


0.554 


WTs’  UTb’ 
MDegK/S  MDegK/S 


0.0088 


-0.0334 
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3.3  ggTnnling  and  Representativeness 


All  of  the  meastirements  needed  to  completely  describe  the  phenomenon  of 
interest  are  rarely  available;  consequently,  the  true-population  statistics  ^e  never 
absolutely  known.  However,  if  a  sufficient  number  of  measurements  are  taken,  it 
is  possible  to  form  acceptable  estimates  of  the  means,  v^ances,  and  other  popula¬ 
tion  statistics.  For  a  statistically  stationary  process,  it  is  possible  to  determine  toe 
number  of  random  independent  samples  needed  to  estimate  ffie  process  popul^on 
statistics  to  within  a  specified  degree  of  precision  and  statisti^  sigmficance.  The 
number  of  samples  required  varies  roughly  with  the  mverse  of  the  square  of  the 
required  precision.  Procedures  for  calculating  the  precision  ^d  signific^ce  are 
available  in  any  number  of  statistics  texts.  Lenschow  et  al.  (1994)  provide  a 
statistical  significance  calculation  procedure  adapted  for  atmospheric 
measurements. 


The  difficulties  in  applying  standard  statistical  procedures  to  meteorological 
time  series  data  are  that  meteorological  processes  are  seldom  stationapr,  and 
sequential  samples  taken  firom  a  time  series  are  neither  r^dom  nor  statistic^y 
independent.  Judgment  on  the  part  of  the  analyst  is  needed  to  overcome  first 
problem.  Means  or  standard  deviations  calculated  firom  a  data  set  obtamed  during 
a  significant  trend  or  transition  may  be  mathematically  correct  but  functionally 
unrepresentative.  For  example,  a  mean  wind  direction  measured  d^ng  a  win 
shift  will  not  represent  the  prevailing  conditions  before  or  after  that  event,  and 
the  standard  deviation  of  these  measurements  wiU  be  excessively  large.  A 
procedure  for  estimating  statistical  independence  is  described  below. 


The  net  effect  of  sequential  sampling  fi*om  a  time  series  where  persistence 
occurs  on  time  scales  greater  than  the  sampling  interval  is  that  the  effective 
(statistically  independent)  sample  size  (n)  is  less  than  the  number  of  data  pomts 
measured  (m).  A  measure  of  persistence  (Ra)  given  by  Brooks  and  Camithers 

(1953) is 


Ra  = 


-1, 


(3-16) 


where  a  is  the  standard  deviation  of  a  data  set  and  Oa  is  the  s  Ward  deviation  of 
the  differences  obtained  from  one  observation  to  the  next  within  the  data  set.  It 
the  data  are  fi-om  a  random  series,  aj  should  approach  a^|2  and  Ra  shoffid 
approach  zero.  An  estimate  of  n  can  be  obtained  by  dividing  m  by  (|Ra|  +1), 
where  |Ral  denotes  the  absolute  value  of  Ra. 
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3.4  Ttj»n1-Time  Predictors 


Predictor  programs  that  forecast  the  next  data  point  in  a  time  sequence 
based  on  past  data  points  can  be  useful  real-time  QC  tools.  Predictors  are  based 
on  the  premise  that  time-sequenced  meteorological  conditions  follow  trends  and 
that  each  data  point  is  correlated  to  its  near  neighbors.  A  combination  of  the 
mean  of  neighboring  data  points  (the  local  mean)  and  ^e  last  valid  measured  data 
point  are  reasonable  estunates  for  the  next  data  point  in  a  time  series.  Predictor 
programs  use  this  estimate  to  create  a  forecast  of  this  data  point.  Bach  new  data 
point  is  then  compared  to  the  predicted  value.  If  the  new  data  point  is  witlm 
boxmds  (within  upper  and  lower  departure  limits  of  the  predicted  value),  it  is 
accepted  without  qualification.  If  the  new  data  point  is  out  of  bounds,  it  is 
flagged.  The  flags  draw  the  attention  of  an  analyst  who  must  then  decide  whether 
each  out-of-boimds  data  point  represents  a  significant  d^ge  in  meteorological 
conditions  or  the  emergence  of  an  instrument  fault  condition.  For  variables  that 
are  expected  to  change  slowly  with  time  (pressure,  for  example),  most  of  the 
weighting  may  be  on  the  running  mean.  For  less  conservative  variables  (wind 
speed,  for  example),  more  weight  may  be  on  the  latest  recorded  wind  speed. 
Weighting  schemes  and  departure  limits  require  "tuning"  for  each  site  and 
measured  variable  to  minimize  the  false-alarm  rate. 

Hojstrup  (1993)  describes  several  real-time  statistical  predictor  models  that 
work  weU  on  data  distributions  that  are  near-Gaussian.  Difficulties  arise  with 
this  technique  when  the  characteristics  of  the  time  series  change.  An  adaptive 
discrimination  factor  can  be  used  at  the  expense  of  additional  computational 
complexity.  This  screening  procedure  is  also  useful  for  ofiT-line  quality  control 
where  the  analyst  has  an  opportunity  to  tune  the  predictor  to  the  data  set. 
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CHAPTER  4 


POST-PROCESSING  QUALITY  CONTROL 


4.1  Post-Processing  Steps 

Post  processing  can  include  a  number  of  analytical  and  statistical 
procedures  that  are  too  computationally  demanding  or  time  consuming  for  use  in 
real-time  QC.  Post-processing  steps  can  include  the  use  of  despikere  to  identify 
out-of-tolerance  data  points,  interpolators  to  fill  in  missing  data  points,  objective 
analysis  of  the  data  field,  spectrum  analysis  to  define  the  contributions  to  variance 
by  different  frequency  intervals,  and  smoothers  that  present  general  trends  with¬ 
out  high  frequency  content.  Post  processing  is  generally  necessary  to  transform 
data  into  information  that  can  be  used  to  model  desired  phenomena.  Each  post¬ 
processing  step  affects  the  information  content  of  a  data  set  and  should  be  done 
with  a  clear  understanding  of  how  the  derived  information  is  to  be  used.  The 
range  customer  should  be  provided  with  a  description  of  the  post-processing  steps 
taken,  the  rationale  for  using  each  procedure,  and  the  likely  effect  on  the  data  set. 


4.2  Despikers 

Noise  can  be  caused  by  power  surges,  radio-fi*equency  interference,  random- 
bit  errors  during  data  transmission,  or  by  loose  or  failing  connections  or  compo¬ 
nents  within  the  sensor  or  transmitter.  These  problems  may  cause  impulse  noise 
or  spikes.  Spikes  are  often  Oaut  not  always)  randomly  distributed,  can  be  of  either 
sign,  and  are  usually  of  short  duration  and  characterized  by  rapid  rise  and  fall 
times. 

4.2.1  Central-Tendency  Measures.  Spikes  are  usually  seen  as  departures  fi-om 
the  central  tendency  of  the  majority  of  points  within  a  data  set.  Therefore,  the 
most  obvious  way  to  identify  spikes  is  to  compare  eadi  data  point  to  a  measure  of 
central  tendency.  Visually  identifiable  spikes  depart  firom  central  tendency  by 
some  user-defined  criterion  such  as  3  standard  deviations  beyond  the  mean.  The 
QC  analyst  should  try  to  imderstand  the  reasons  spikes  may  be  present  and  how 
they  are  distributed  within  the  data  before  trying  to  remove  those  spikes. 

Spikes  that  are  clustered  or  grouped  present  a  much  more  difficult  QC 
problem  than  those  that  are  randomly  ^stributed  through  a  data  set.  Filter 
lengths  must  be  carefully  chosen  to  exceed  the  length  of  any  spike  cluster. 
Despikers  work  reasonably  well  when  spikes  are  clearly  distinguishable  fi-om  valid 
readings,  when  the  number  of  spikes  is  10  percent  or  less  of  the  total  data  set,  and 


4-1 


when  the  spikes  are  not  clustered.  Clustered  spikes  present  a  spea^  problem 
because  an  adaptive  filter  can  adjust  to  the  spike  cluster  and  bepn  flag^g  valid 

data  points.  Discrimination  between  an  abrupt  sl^  ^ 

emergence  of  clustered  spikes  is  one  of  the  most  diffic^t  QC  problems.  Some 
degree  of  mannal  intervention  by  the  analyst  will  likely  be  required. 

The  available  measures  of  central  tendency  are  the  mean,  median,  and 
mode  The  mean,  or  average,  is  the  most  famibar  and  the  most  firequently  used 
for  data  processing,  because  it  is  easy  to  compute.  The  median,  representmg  the 
50th  percentile  point  in  a  distribution,  is  less  common  because  it  requires  a 
sorting  procedure  that  is  more  computationally  intensive  ^an  computetion  of  the 
mean  The  mode,  the  most  fi-equently  observed  data  pomt  value  m  the  ^ta  set,  is 
a  useful  measure  of  central  tendency  if  the  magnitude  of  the  ^e  cental  tendency 
(or  conversely,  the  magnitude  of  the  spike)  is  constant  oyer  the  sainphng  period. 
Mode-based  arguments  are  used  to  define  "consensus  with  radar  wmd  profiler 

data. 

4.2.2  Mean  Filters.  A  variety  of  mean  filters  are  available  for  identifying  spkes 
in  a  data  set.  The  simplest  of  these  filters,  described  m  chapter  3,  para^aph  3.2, 
uses  the  data  set  maximum,  minimum,  mean,  and  stand^d  deyiatiom  A  ratio  ot 
range  to  standard  deviation  is  computed.  If  this  ratio  is  l^ge,  u^catmg  the 
probable  presence  of  spikes,  a  criterion  such  as  3  standard  deviations  beyond  the 
mean  is  used  to  flag  the  spikes.  While  this  simple  procedure  is  useM  for  means 
that  remain  relatively  constant,  the  means  in  meteorological  data  often  vary 
greatly  with  time.  A  data  point  in  a  time  series  usually  exhibits  a  d^er  relation¬ 
ship  to  its  near  neighbors  than  to  other  points  that  are  more  dis^t  fi*^  it. 
(Recurring  periodic  oscillations  are  exceptions  to  this  generahzation.)  Therefore, 
local  or  running  mean  is  usually  a  more  effective  measure  of  central  tendency  for 
spike  detection  in  time-series  data. 

Roberts  (1993)  describes  a  running-mean  filter  used  by  NASA  to  filter 
spurious  points  from  time  series.  A  running  mean  of  100  data  points  is  estab- 
fished  at  the  beginning  of  each  data  set  with  each  subsequent  dat^oint  <»mpared 
to  the  mean  and  standard  deviation  of  the  preceding  100  pomts.  (These  100 
points  are  first  examined  to  ensure  that  they  are  relatively  fi-ee  of  spikes  )  the 
data  point  falls  within  the  user-defined  criterion  (usually  within  4  or  5  st^dard 
deviations  of  the  mean),  the  point  is  accepted  into  the  mean.  The  filter  then  steps 
forward  one  increment  in  time,  dropping  its  earliest  data  pomt,  adding  its  newly 
accepted  data  point,  and  beginning  the  process  again  by  rewmputmg  a  new  mean 
and  standard  deviation.  Roberts  (1993)  includes  a  FORTRAN  program  for  this 
running-mean  filter. 

Running-mean  filters  of  the  Roberts  (1993)  design  offer  an  effectwe 
methodology  for  spike  identification.  It  can  run  through  single  or  multiple  passes. 
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as  necessary,  to  clean  up  a  noisy  data  set.  One  disadvantage  of  this  technique  is 
that  it  is  fairly  computer  intensive,  requiring  computation  of  both  the  mean  and 
standard  deviation  at  each  step.  Another  disadvantage  for  meteorological  date 
applications  is  that  a  large  number  of  data  points  is  needed  in  the  filter  (on  the 
order  of  100)  to  compute  a  statistically  stable  standard  deviation.  A  simplification 
of  this  technique  would  be  to  use  a  standard  deviation  representative  of  the  whole 
date  set  and  recompute  only  the  means.  This  tec^que  will  work  so  long  as  the 
standard  deviation  remains  reasonably  invariant  in  time.  Alternative  running- 
mean  techniques  include  adaptive-running  mean  with  threshold  logic  and  clipped- 
mean  filters. 

An  adaptive  running-mean  filter  with  threshold  logic  and  dual  levels  of 
flagging  was  developed  for  meteorological  date  applications  by  Biltoft  (1993).  The 
adaptive  filter  adjusts  to  abrupt  changes  in  signal  strength  using  dual-level 
flagging  threshold  logic  and  a  filter  reset  designed  for  adjustment  to  abrupt 
changes  in  signal  level.  Flagging  criteria  are  based  on  user-chosen  levels  of 
departure  fi-om  the  running  mean.  These  departure  levels  are  usuafiy  based  on  a 
standard  deviation  estimate.  One  level  of  flagging  is  chosen  for  easily  identified 
spikes  such  as  date  points  that  are  4  or  5  standard  deviations  fi-om  the  running 
mean.  These  unconditionally  flagged  points  are  excluded  fi-om  the  running-filter 
mean.  The  second  flagging  level,  set  within  the  range  of  the  first  (2  or  3  standard 
deviations,  for  example),  conditionally  flags  a  suspect  datum.  This  point  is 
included  within  the  running-mean  filter  and  is  re-examined  upon  exiting  the  rear 
of  the  filter.  If,  in  the  meantime,  the  filter  has  adjusted  to  within  the  accepted 
range  of  the  datum,  the  flag  is  removed  and  the  datum  is  wnsidered  valid.  If  the 
filter  has  not  adjusted  to  within  range,  the  datum  retains  its  flag.  Because  it  does 
not  require  continuous  standard  deviation  computations,  it  is  relatively  effident 
computationally  and  can  use  short  running-mean  filter  lengths  (10  to  25  points) 
which  are  more  suitable  for  highly  variable  date  sets  with  short  autocorrelation 
times.  This  technique  works  reasonably  well  on  noisy  date  sets  as  long  as  spike 
clusters  do  not  exceed  filter-reset  length.  Results  of  the  application  of  this 
filtering  technique  to  the  date  presented  in  figure  3-1  is  shown  in  figure  4-1. 


Clipped-mean  or  trimmed-mean  filters  are  running-mean  filters  formed 
from  a  number  of  date  points  on  either  side  of  the  datum  in  question.  A  user- 
selected  number  (usually  1  or  2)  of  the  highest  and  lowest  values  in  the  set  are 
removed  (clipped),  and  a  new  mean  is  computed.  If  the  datum  in  question  is 
among  the  points  removed,  it  is  flagged  or  replaced  by  this  mean.  If  not,  the 
datum  is  accepted  as  valid.  This  technique  can  be  computationally  more  efficient 
than  the  previously  mentioned  running-mean  filter,  but  it  is  more  susceptible  to 
spike  clusters.  As  the  munber  of  clipped  points  is  increased,  the  clipped  mean 
approaches  the  median  as  a  limit.  Restrepo  and  Bovik  (1988)  describe 
applications  of  an  adaptive  trimmed-mean  filter  for  reducing  noise  in  digital 
signals. 
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Figure  4-1.  Edited  v-component  sonic  anemometer  time  series  data  from  instruments  moimted  at  the  2-, 
4-,  6-,  and  8-meter  levels  on  a  tower. 


4  2  3  Median  Filters.  Median  filters  provide  an  effective  means  of  separatmg 
signal  fi-om  noise  provided  that  the  filter  is  large  enough  to  ^scnminate  grouped 
spikes.  Median  filters  operate  by  comparing  the  observed  value  m  the  ^nter  of  a 
sliding  window  with  the  median  of  data  points  that  form  the  wmdow.  Window 
length  (W)  of  a  median  filter  is  determined  by  its  order  (the  number  of  data 
points,  N,  on  each  side  of  the  central  point)  defined  by  the  relation  (Gallager  and 

Wise,  1981). 


W  =  2N  +  1 


(4-1) 


A  first-order  (N=l)  filter  has  a  window  length  of  3  (the  center  pomt  ^  data 
points  on  each  side);  a  second-order  filter  has  a  window  length  of  5.  'me  choice  of 
N  for  a  median  filter  is  crucial.  It  cannot  discriminate  against  any  pulse  or  pike 
cluster  whose  width  is  greater  than  N.  Median  filters  are  particularly  pefiil  with 
very  noisy  data  sets  as  long  as  N  is  sufficiently  large  to  discriminate  spike 
clusters.  A  running  median  can  be  less  sensitive  to  noise  than  a  running  mean  of 
the  same  order,  because  it  is  influenced  by  the  number  but  not  by  the  magmtude 
of  noise  spikes.  However,  if  the  number  of  spikes  sampled  exceeds  N,  the  medip 
will  represent  the  spikes  rather  than  the  noise-fi:ee  data.  A  median  filter  should 
track  either  the  valid  data  or  the  spikes  but  (unlike  a  mean  filter)  not  some 
average  of  the  two. 

Use  of  threshold  logic  with  a  median  filter,  as  described  by  Brock  (1986), 
requires  two  passes  through  the  data  set.  'The  results  of  the  first  pass  are  used  to 
construct  a  histogram  fi-om  the  difference  values  between  the  data  and  ffiter 
outputs.  Spikes  appear  as  lobes  in  what  would  otherwise  be  a  normal  distribu¬ 
tion.  The  threshold  is  defined  as  the  minimum  on  the  histogram  that  separates 
the  normal  distribution  of  "good"  data  points  firom  lobes  caused  by  spikes.  The 
second  pass  through  the  data  set  targets  these  lobes.  Additional  passes  can  be 
made,  as  necessary,  to  converge  on  a  smoothed  signal. 

The  primary  disadvantage  of  median  filters  is  the  computation^y 
demanding  sorting  needed  to  define  the  median.  As  computer  power  increases, 
computational  requirements  become  less  of  an  obstacle.  Astola  and  Cmnpbell 
(1989)  present  a  fast  median-filter  computation-algorithm  based  on  a  "double 
heap"  process. 

4.3  Interpolation  Procedures 

Interpolation  is  used  to  replace  bad  or  missing  data  points  with  estimates 
that  closely  reflect  the  real  conditions  at  the  time  of  measurement.  Interpolation 
is  not  always  desirable.  It  is  not  necessary  if  data  are  to  be  used  simply  to 
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generate  summary  statistics.  However,  inteipoktion  may  be  needed  ^ 

fet  is  to  be  subjected  to  spectrum  analysis  or 

time  series  with  no  gaps.  Sometimes  a  "best  guess  T 

information  gap.  The  customer  should  be  consulted  concermng  the  desirabihty  of 

using  interpolated  data. 

Simple  interpolation  for  one  or  two  missing  sequential  (^ta  pomts  can  be 
done  using  the  average  of  the  readings  before  and  after  the  mssmg  datnm  T^s 
simple  procedure  is  based  on  the  observed  autocorrelation  withm  .“^eteorologic^ 
data  Each  data  point  is  more  closely  related  to  its  immediate  neighbors  th^  to 
points  at  a  greater  distance  in  the  time  series.  The  rmmmg  me^  and  median 
mters  described  in  paragraph  4.2  can  also  be  used  to  fiU  m  mssmg  ‘iate  pomt^ 
although  these  filters  produce  interpolated  values  closer  to  the  cental  ^^^en^ 
Bednar  (1983)  describes  a  procedure  where  the  datim  m^estion  is  replaced  by 
the  median  of  a  sliding  window  of  near-neighbor  pomts.  The  mean  of  » 
unaltered  if  it  is  used  to  fill  in  missing  data  points,  but  this  replacement  is  done 
at  the  expense  of  local  data  point  autocorrelation. 


4.4  Ohiective  Analysis. 

Objective  analysis  is  a  procedure  used  to  perform  QC  on  date  fields 
consisting  of  multiple  points  in  an  array.  Each  datum  is  comp^ed  te  its  ne^ 
neighbors.  This  procedure  is  particularly  useful  for  networ^  of  surface  stetions 
where  known  physical  relationships  (the  hypsometric  equation,  for  exampl^  can 
be  used  to  form  estimates  against  which  measured  values  are  comp^ed.  The  end 
product  is  a  smoothed  distribution  of  values  that  are  mutuaUy  consistent  with 
each  other  and  with  known  physical  relationships.  Wade  (1987)  describes  an 
objective  analysis  technique  used  on  surface  mesomet  pressure,  temperature, 
humidity,  and  wind  date. 


4.5  Spectrum  Analysis 

Spectrum  analysis  is  perhaps  the  most  powerful  tool  available  for  quality 
control.  Fault  conditions  that  are  transparent  in  the  time  doinain  ^e  often  de^ly 
revealed  in  the  fi-equency  domain.  UntU  recently,  Fourier  analysis  has  donated 
frequency  domain  applications;  however,  wavelet  analysis  is  now  also  available  as 

a  quality  control  tool. 


4  5  1  Fourier  Analysis.  Measurements  made  at  regular  intervals  are  often  pre¬ 
sented  as  time  series.  Time-series  plots  such  as  figure  3-1  and  4-1  representing 
crossaxis  (v-component)  wind  measurements  fi-om  four  some  ^emometers  permit 
a  visual  scan  of  date  fields  for  interesting  features  and  possible  faults  such  as 
missing  data  or  spikes.  Statistical  analyses  of  time  series  also  provide  information 
about  the  mean,  variance,  and  covariance.  Time  series  analysis  alone  cannot 
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describe  how  the  variance  of  a  tune-var3dng  quantity  is  distributed  among  the 
various  scales  of  motion  that,  when  combined,  form  the  time  series.  Xo  better 
imderstand  atmospheric  processes  and  to  eliminate  unwanted  or  erroneous 
signals,  it  is  often  necessary  to  transform  the  series  from  the  tune  domain  to  the 
frequency  domain  where  spectral  energy  is  partitioned  by  frequency  rather  than 
by  time.  Fourier  analysis  is  based  on  the  premise  that  an  extended  time  series 
pan  be  decomposed  into  a  linear  combination  of  exponentials  that  are  represen¬ 
table  as  a  spectrum  of  sinusoidal  components.  The  basic  Fourier-based  power- 
spectrum  analyses  are  presented  by  Blackman  and  Tukey  (1958).  For  many  kinds 
of  geophysical  phenomena  (the  flow  of  tides,  daily  or  annual  temperature  change), 
Fourier  analysis  cjan  be  used  to  define  the  most  energetic  and,  therefore,  the  most 
important  frequencies  in  a  time  series  and  to  illustrate  details  of  the  interactions 
between  these  phenomena. 

Fourier  analysis  of  time-series  data  can  provide  statistical  measures  of 
instrument  performancje  and  insight  into  the  physics  governing  the  interactions 
between  sets  of  time  series  data.  Simple  linear  correlation  and  multivariate- 
analysis  procedures  are  not  appropriate  for  time-series  data,  because  significance 
tests  are  based  on  the  hypothesis  that  the  variate  samples  are  imcorrelated. 
Autocorrelation  is  nearly  always  present  in  the  time  series  of  geophysical  data. 
Also,  the  relationship  between  two  time  series  can  be  obscured  by  an  in-phase 
relationship  at  some  frequencies  and  an  out-of-phase  or  phase-lagged  relationship 
at  other  frequencies.  Multivariate  analyses  can  be  applied  to  the  spectra  of  a  time 
series  after  a  Fourier  transformation  of  the  series  from  the  time  domain  to  the 
frequency  domain.  Transformation  partitions  the  power  (energy  expenditure  per 
unit  time)  of  the  time  series  into  harmonic  frequency  components.  For  a  given 
variable,  "power"  is  equivalent  to  the  change  in  that  variable  per  umt  time.  The 
total  power  of  the  process  is  equal  to  the  sum  of  the  contributions  by  the  harmonic 
components.  The  advantage  of  power-spectrum  analysis  over  the  usual  multi¬ 
variate  analysis  is  that  the  amoimt  of  power  contributed  by  one  harmonic  is 
independent  of  the  amplitudes,  phases,  and  frequencies  of  the  other  harmonics  in 
the  time  series.  Detailed  information  on  Fourier  transform  statistical  applications 
is  foimd  in  texts  on  time-series  analysis  (Koopmans,  1974;  Jenkins  and  Watts, 
1968);  general  information  on  time-series  analysis  for  meteorological  data  are 
presented  by  Stull  (1988). 

Relationships  between  two  time  series  can  be  evaluated  by  examination  of 
the  spectra  from  the  two  series  and  their  complex  products,  the  cross  spectrum, 
phase  (PHASE),  and  coherence  (COH).  The  cross  spectrum  is  represented  by  the 
cospectrum  (COS)  and  quadrature  (QUAD).  The  cospectrum  is  analogous  to  an  in- 
phase  covariance  between  the  two  spectral  frequencies,  and  the  quadrature 
spectrum  is  a  similar  measure  phase-shifted  wavelength  (90°).  The  COS  and 
QUAD  define  the  cross-spectrum  covariance  components  in  a  Cartesian-like 
coordinate  system,  while  PHASE  and  COH  express  these  components  in  a 
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nonnEdized  polar  coordinate  system.  The  PHASE,  a  representation  of  the  angular 
relationship  between  two  spectra,  is  defined  as 


PHASE  =  ATAN2{QUADICOS), 


(4-2) 


where  ATAN2  is  defined  in  FORTRAN  as  the  arctangent  fimction  expanded  to  the 
range  -%  to  n.  The  PHASE  is  presented  in  degrees.  The  COH,  the  squared 
coefficient  of  coherence,  is  a  measure  of  the  correlation  between  two  time  senes  as 
a  function  of  firequency  given  by 


COH  =  (COS^  +  QUAD^KSPCTR1)*(SPCTR2), 


where  SPCTRl  and  SPCTR2  are  the  spectral  component  magnitudes  of  time  series 
1  and  2.  The  COH  is  dimensionless  and  ranges  in  magnitude  fi*om  0  to  1. 
Statistical  significance  tests  can  be  applied  to  coherence  data  using  procedures 
presented  in  Koopmans  (1974). 

A  Fast-Fourier  transformation  (FFT)  and  spectrum  analysis  program 
designed  for  use  with  meteorological  data  is  described  by  Kaimal  and  Gaynor 
(1983).  The  spectra  produced  by  this  FFT  are  smoothed  through  block  averagmg, 
tapered,  scaled  to  meter-kilogram-second  (mks)  units,  and  normahzed  by  multi¬ 
plying  each  harmonic  component  by  its  fi*equency.  For  example,  table  4.1  presents 
alongaxis  (u-component)  wind  spectra  obtained  fi’om  sonic  anemometers  moimted 
at  0.5-  and  1.0-meters  above  ground  level  on  a  tower.  Measurements  fi'om  the 
sonics  were  taken  concurrently  with  SPCTRl  representing  the  0.5-meter  level  and 
SPCTR2  representing  the  1.0-meter  level.  Spectrum  band  centroid  fi-equencies  are 
presented  in  the  left  column,  and  coherence  between  the  spectrum  fi*equency 
bands  is  shown  in  the  right  column.  Coherence  between  the  0.5-  and  1.0-meter 
components  is  high  at  the  lower  fi*equencies  (0.019  to  0.063  Hz)  but  dimimshes 
rapidly  at  the  higher  frequencies.  Analyses  of  this  type  can  be  used  for  quahty 
control  to  locate,  by  frequency,  noise  components  in  a  data  set. 

4  5  2  Wavelet  Annlvsis.  Fourier  analysis  is  based  on  the  notion  that  an  infinite 
time  series  can  be  represented  in  the  frequency  domain  as  a  spectrum  of 
sinusoidal  components  of  varying  frequencies.  However,  some  geophysic^ 
phenomena  are  characterized  by  intermittent  bursts  of  energy  or  discontinuities 
rather  than  by  regular  oscillatory  motions.  Turbulence,  for  example,  occurs  as  a 
series  of  discontinuities  in  an  otherwise  quiescent  background.  On  a  larger  scale, 
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TABLE  4-1.  SONIC  TEMPERATUEE  TIMES  SERIES  SPECTRA,  COSPECTRA,  (QUADRATURE,  AND  THE 
PHASE  AND  COHERENCE  BETWEEN  THESE  TIME  SEEJES. 

SPECTRA  DATA  SU3VIMARY 
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abrupt  changes  in  the  wind,  pressiire,  and  temperature  fields  associated 
frontal  passage  are  discontinuities  that  disrupt  an  otherwise  regular  oscillating 
pattern.  Fourier  analysis  is  not  weU  suited  for  investigation  of  these  intermittent 
processes,  because  it  does  not  localize  the  event  in  time,  and  it  requires  multiple 
frequencies  to  represent  discontinuities  or  intermittencies. 

The  limitations  of  Fourier  analysis  have  led  to  the  development  of  a 
technique  known  as  wavelet  analysis.  Wavelets  are  zero-mean  functions 
consisting  of  short  oscillations  locahzed  in  both  time  and  frequency.  A  wavelet  is 
characterized  by  a  dilation  factor  which  affects  its  size  and  amplitude  and  by  a 
translation  parameter  which  defines  its  origin  or  position  within  a  time  senes. 

The  wavelet  transform  is  a  convolution  (the  inverse  transform  of  a  cross-spectrum) 
of  a  wavelet  with  a  time  series.  The  resultant  covmance  product  localizes  small- 
scale  features  with  fine  spatial  resolution,  permitting  a  detailed  analysis  of 
intermittent  phenomena  that  is  not  possible  with  traditional  Fourier  analysis. 

Wavelet  analysis  is  useful  in  locating  intermittencies  and  discontinuities, 
because  the  product  of  the  wavelet  with  time-series  elemente  is  greatest  where  ^e 
change  in  time-series  element  magnitude  is  greatest.  The  time  series  illustrated 
in  figures  3.1  and  4.1  were  subject  to  a  discrete  wavelet  transform  with  the 
Lemarie-Meyer-Battle  (LMB)  wavelet  basis  using  procedures  developed  by 
Kosteniuk  (1993).  Inverse  wavelet  transforms  were  then  performed  to  separate 
the  signal  content  by  wavelet  scale  indices.  Figures  4.2  and  4.3  present  the 
imedited  and  edited  v-component  signal  components  at  scale  index  13 
(corresponding  roughly  to  a  frequency  of  2.5  Hz).  Figure  4.2  contains  large 
signatures  for  the  data  at  4-,  6-,  and  8-meter  levels  correspon^g  to  the  noise 
spikes.  The  much-reduced  signature  amplitudes  in  figure  4.3  indicate  that  the 
large  spikes  have  been  removed  by  the  despiking  process. 

The  field  of  wavelet  analysis  has  essentially  developed  since  the  1980s  with 
primary  applications  in  optics,  acoustics,  and  signal  compression.  A  nmthematical 
description  of  wavelets  is  presented  by  Meyer  (1992).  Among  geophysical 
applications,  wavelet  analysis  is  being  used  to  describe  and  model  turbulence 
intermittencies  (Farge,  1992;  CoUineau  and  Brunet,  1993).  Because  it  offers  the 
promise  of  locahzing  signal  components  by  frequency  and  time,  wavelet  analysis 
should  evolve  into  a  powerful  quality  control  tool. 
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Figiire  4-2.  Unedited  v-component  sonic  anemometer  time  series  at  wavelet  scale  index  13.  (Note:  The 
entire  figure  is  autoscaled.) 
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4.6  Smoothers 


In  some  situations,  a  smoothed  average  value  or  trend  is  more  important 
than  point-to-point  detail  which  may  contain  noise  or  other  undesirable 
characteristics.  Smoothers  remove  unwanted  high-jfrequency  information  to 
expose  averages  and  trends  that  may  otherwise  be  indiscernible  in  a  fluctuating 
data  field.  Smoothers  operate  on  the  principle  that  noise  is  randomly  distributed 
and  can  be  removed  fi-om  a  data  set  by  simple  averaging.  This  assumption  may 
not  always  be  good,  so  a  data  set  should  be  examined  carefully  before  smoothing. 
Despiking  may  be  needed  prior  to  the  smoothing  operation. 

Smoothing  can  be  achieved  by  simple  block  averaging  or  by  point-by-point 
replacement  of  each  datum  by  the  mean  or  median  of  its  neighbors.  Block 
averaging  compresses  the  data  set  but  may  introduce  discontinuities  or  jumps 
between  successive  nonoverlapping  blocks.  Replacement  of  each  datum  by  a 
measure  of  central  tendency  produces  a  much  smoother  data  set.  Smoothing  is 
used  extensively  for  procedures  like  image  restoration  (Restrepo  and  Bovik,  1988) 
and  for  filtering  out  impulse  noise  (Palmieri  and  Boncelet,  1989).  Smoothing  of 
meteorological  data  is  most  useful  for  conservative  variables  such  as  pressure  or 
potential  temperature  where  consistency  with  known  physical  relationships  is 
more  important  than  the  details  of  high-fi:equency  variations. 
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APPENDIX  A 


INSTRUMENT  FAULT  CONDITIONS 


A1  Mecbanicfll  Anemometers 

Mechanical  anemometers  are  vulnerable  to  a  variety  of  mediamcal  and 
electrical  faults.  These  faults  can  usually  be  identified  by  spot  checking  the  data. 

A.  1.1  Threshold  Errors.  Rotating  mechanical  wind  instruments  are  mounted  on 
low  fiiction  bearings.  These  bearings  eventually  wear  down  and  accumulate  dust 
causing  bearing  fiiction  to  increase.  Because  an  increase  in  bearing  fiiction 
requires  an  increase  in  turning  torque  to  cause  movement,  the  threshold  speed  is 
increased. 

An  increase  in  bearing  fiiction  is  usually  an  imperceptibly  slow  process  that 
is  difficult  to  detect  operationally.  Performing  threshold  speed  checks  in  a  wind 
tunnel  is  the  best  way  to  define  bearing  wear.  A  simple  (but  nonquantitative)  test 
of  torque  described  in  the  draft  Environmental  Protection  Agency  (EPA  QA)  hand¬ 
book  (Lockhart,  1989A)  is  to  roll  the  sensor  housing  slowly  on  a  smooth,  level,  flat 
surface  and  observe  the  shaft.  If  the  shaft  turns  with  the  sensor  housing,  the 
torque  may  not  meet  specification.  A  quantitative  measure  of  torque  on  the 
turning  shaft  can  be  obtained  using  the  R.  M.  Yoimg  torque  disk.  ^  For  field  testing 
of  mounted  sensors,  starting  torque  can  be  measured  using  precision  torque 
instruments  that  fit  onto  the  turning  shaft.  Available  torque  instruments  of  this 
type  include  the  Waters  torque  watch  and  the  Gm-cm  torquemeter.  Sensor 
manufacturers  should  be  consulted  about  the  torque  sensor  that  works  best  with 
their  instruments.  A  crude  but  effective  bearing  chedc  for  an  installed  propeller 
anemometer  is  to  hang  a  small  paper  dip  on  one  blade  and  observe  the  rotation  of 
that  blade  towards  the  groxmd. 

Threshold  effects  are  difficult  to  observe  unless  wind  speeds  are  bracketing 
the  cup-breakaway  speed.  Near  neighbor  and  velodty  profile  checks  can  aid  in 
verifying  anemometer  performance,  but  bearing-drag  problems  can  be  most  clearly 
observed  by  examining  wind-speed  histograms  or  distribution  tables.  Bearing 
drag  should  load  the  "calm"  wind-speed  bin  at  the  expense  of  the  next  higher  bin, 
as  illustrated  by  the  histograms  in  figure  A.I. 

A.  1,2  MpohflniPfll  Anemometer  Overspeeding.  Mechanical  anemometers  require 
a  certain  amoxmt  of  air  to  pass  through  them  before  they  approach  equilibrium 
with  a  new  wind  speed.  This  length  of  air  passage  remains  fairly  invariant  over 
the  sensor’s  linear  operational  range  and  is  known  as  a  distance  constant 
(Lockhart,  1987).  Because  of  their  asymmetric  drag,  cup  anemometers  accelerate 
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faster  than  they  decelerate,  the  result  is  a  net  overspeed  condition  in  turbxilent 
winds.  The  overspeed  ranges  from  1  to  10  percent  depending  on  the  scale  ^d 
intensity  of  turbulence.  Corrections  are  difficult  to  calculate,  so  this  effect  is 
usually  overlooked;  however,  overspeeding  can  be  minimized  by  proper  placement 
of  instrumentation.  A  site-selection  criterion  for  cup  anemometers  should  include 
mininmiTn  exposure  to  turbulent  conditions  such  as  locations  in  the  wake  of  trees 
and  buildings  or  in  the  proximity  of  rough  surfaces.  Helicoid-propeller  anemo¬ 
meters  are  virtually  symmetric  and  exhibit  less  overspeeding  in  turbulent 
conditions.  Propeller  anemometers  may,  therefore,  be  more  suitable  than  cup 
anemometers  for  use  in  areas  subjected  to  strong  fluctuating  wind 
conditions. 

A.  1.3  Sensor  DgmagPi-  The  aerodynamic  performance  of  rotating  anemometers  is 
critically  dependent  on  instrument  balance  and  symmetry.  If  a  blade  or  cup  is 
misshapen,  broken,  or  improperly  mstalled,  the  sensor  will  still  report  data  but 
performance  will  be  degraded.  This  degradation  results  in  uncharacteristic  signal 
fluctuations,  particularly  at  low  wind  speeds.  Sensor  damage  can  be  detected  by 
observing  the  variance  or  standard  deviation  of  wind  speed  readings. 

A.1,4  Battery  Voltage.  A  12-volt  battery  is  often  used  to  power  mechanical 
anemometers  and  vanes.  As  the  battery  discharges  and  voltage  drops,  the 
operating  span  of  the  instruments  decreases  causing  erroneous  data.  When 
battery  voltage  is  used  as  a  reference,  tbis  voltage  should  be  monitored  and 
reported  on  an  additional  data  chaimel.  When  feasible,  a  battery  charger  should 
also  be  stationed  at  the  measurement  site. 

A.  1.5  Negative  Wind  Sneed.  A  rotating  anemometer  may,  tmder  near-calm 
conditions,  report  a  small  negative  wind  speed  (-0.1  m/s,  for  example).  Tffis 
negative  wind  speed  may  be  an  artifact  of  extension  of  the  linear  cahbration  curve 
into  the  sensor’s  nonlinear  region.  In  this  case,  the  instrument  may  not  be  in  a 
true-fault  condition  unless  there  is  supporting  evidence  to  suggest  otherwise. 

A.  1.6  Signal  Check.  An  end-to-end  signal  dbeck  for  a  rotating  anemometer  can 
be  accomplished  in  field  conditions  by  removing  the  sensor  head  and  installing  a 
small  motor  driven  shaft  that  rotates  at  a  known  rate.  The  output  signal  at  the 
data-recording  device  is  compared  to  the  output  expected  for  that  rotation  rate. 
This  signal  check  verifies  the  performance  of  everything  except  the  sensor  head. 

A.2  Wind  Vanes 


Electrical  or  mechanical  faults  that  occur  with  wind  vanes  are  usually 
associated  with  environmental  degradation  or  animal-induced  damage.  The  tail 
assembly  is  the  component  most  susceptible  to  damage.  Metal  tail  assemblies. 
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which  offer  greater  damage  resistance  at  the  expense  of  responsiveness,  may 
needed  at  locations  where  tail  damage  is  likely  or  at  infrequently  visited  remote 

sites. 

A  2  1  Potentiometer  Errors.  Dust  or  corrosion  can  cause  dead  spots  or  shorts  on 
potentiometer-based  wind  equipment.  A  field  check  for  this  condition  involves 
slowly  rotating  the  vane  both  clockwise  and  counterclockwise  through  its  entire 
electrical  span  while  recording  output.  Note  that  wind  vanes  with  two  potenfio- 
meters  have  a  540°  span.  Changes  in  wind  direction  output  should  be  smooth  and 

continuous.  A  dead  spot  appears  as  a  discontinuity  in  the  reading.  For  plotted 

wind-vane  data,  a  dead  spot  may  appear  as  an  unusual  accumulation  of  data 
points  within  a  very  narrow  range  as  the  wind  vane  moves  through  the  affected 

wind-direction  section. 

Wind-vane  performance  is  not  well  characterized  by  means  and  sl^dard 
deviations,  particularly  for  data  sets  collected  over  extended  periods  of  time. 
However,  information  can  be  gained  by  examining  histogr^s  or 
number  of  wind  readings  per  wind-direction  interval.  A  histogram  with  36  bins 
(10°  intervals)  will  usually  be  sufficient  to  identify  faults  with  a  vane  s  penor- 
mance  or  exposure.  The  analyst  should  examine  bins  that  contain  unusually  high 
or  low  numbers  of  wind  readings.  These  wind-reading  counts  could  mdicate  a 
problem  with  the  equipment,  with  siting  or  exposure,  or  with  the  occurrence  ot 
interesting  terrain-induced  wind-flow  patterns. 

A  2  2  Bent  or  Dam  a  gad  Vane.  Wind  vanes  are  susceptible  to  mechamcal  damage 
caused  by  environmental  factors  such  as  birds  perching  on  the  arms  and  of  ™d 
and  radiation  damage  to  the  aerodynamic  tail.  A  bent-v^e  arm  can  cause  the  tail 
to  hang  up  on  its  supporting  mount,  which  could  appear  in  a  histogram  plot  as  an 
unusual  absence  of  readings  from  certain  wind  directions  coupled  with  an  accumu¬ 
lation  of  readings  from  the  point  where  the  vane  motion  restriction  occurs. 


Damage  to  the  aerodynamic  tail  of  a  vane  may  cause  subtle  changes  in  v^e 
threshold,  distance  constant,  and  damping  ratio  that  are  unlikely  to  be  detected  in 
the  field  Bent  tails  will  impart  a  measurement  bias  and  should  be  replaced,  it  a 
vane  is  used  to  generate  wind  variance  information  and  noticeable  portions  of  the 
tail  are  broken,  the  vane  should  be  replaced.  Wind  vanes  used  to  obtain  mean- 
direction  information  only  can  perform  reasonably  well  unless  a  major  portion  of 
the  tail  is  missing. 

A  2  3  Vgr>P  Tmhalflnrfi-  Wind  vanes  are  balanced  to  align  the  balance  point  over 
the  sensor’s  center  shaft.  Gross  vane  imbalance  can  cause  the  vane’s  toshold 
and  damping  ratio  to  deviate  from  manufacturer’s  specifications,  resulting  m 
unusually  high  or  low  wind  direction  variances. 
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A.3  Thermometers 


Measurement  representativeness,  rather  than  mstrument  malfunction,  is 
the  major  problem  with  thermometers.  Thermometer  sensors  should  be  isolated 
as  much  as  possible  from  electrical  and  environmental  interferants. 

A.3.1  Temperature  Measurement  Problem.  Most  temperature  measurement 
sensors  used  for  meteorological  applications  are  of  the  immersion  type,  which 
includes  probes  that  absorb  energy  from  or  emit  energy  to  their  surroundings. 
These  thermometers  are  calibrated  in  stable  thermal  baths  or  chambers 
maintained  by  laboratories  that  can  reference  these  baths  to  an  acceptable 
temperature  standard.  Once  cahbrated,  a  thermometer  will  likely  remain 
accurate  imtil  external  influences  such  as  foreign-subst^ce  deposits  or  corrosion 
compromise  the  sensor  housing,  or  thermal  and  mechamcal  stresses  alter  the 
thermometer  sensor.  As  long  as  a  thermometer  probe  is  not  subject  to  unrea¬ 
sonable  stress  and  the  system  is  calibrated  at  suitable  intervals,  accuracy  is 
usually  not  a  problem. 

The  major  problem  with  atmospheric-temperature  measurements  is 
representativeness.  A  condition  in  which  the  probe  is  in  thermal  equihbrium  with 
the  surrounding  air  is  desired  for  meteorological  temperature  measurements. 
However,  the  air  is  fairly  transparent  to  radiation  and  changes  temperature 
primarily  through  convection,  while  the  thermometer  probe  surface  responds  to 
solar  radiation,  convection,  and  conduction.  Because  the  absorptance,  conduc¬ 
tance,  aTiH  emittance  of  the  probe  and  the  air  are  different,  the  air  and  probe 
respond  differently  even  when  exposed  to  the  same  thermal  and  radiative  envi¬ 
ronment.  The  result  is  that  a  thermometer  probe  mounted  in  the  atmosphere  is 
rarely  in  equilibrium  with  the  surrounding  air.  Consequently,  a  thermometer  can 
meet  all  criteria  for  accuracy  and  still  provide  measurements  that  are  not 
representative  of  the  true  air  temperature.  Luers  (1990)  describes  a  methodology 
for  estimating  the  temperatxzre  error  on  a  common  temperatxire  probe,  the  rod 
thermistor. 

A.3.2  Radiation  Effects.  Thermometer  probes  absorb  and  emit  radiation  at  a 
much  greater  rate  than  the  surrounding  air.  A  typical  unshielded  probe  will  read 
temperatures  higher  than  the  surroxmding  air  when  exposed  to  direct  sunlight 
during  the  day  and  lower  temperatures  at  night.  Probes  are  often  installed  in 
radiation  shields  to  minimize  this  effect.  Radiation  shields  with  a  ventilation  rate 
of  3  to  5  m/s  can  reduce  radiation  errors  to  several  tenths  of  a  degree.  Radiation 
shields  should,  if  possible,  be  aspirated  mechanically.  Nonmechanical  (naturally 
aspirated)  shields  are  sometimes  used  at  remote  sites  where  power  availability  is 
limited.  Measxirement  representativeness  is  degraded  for  sensors  mounted  in 
naturally  aspirated  shields  dtiring  light  wind  conditions.  Figure  10  in  Wade 
(1987)  illustrates  the  magnitude  of  the  temperature  error  that  can  exist  for  a 
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typical  temperature  probe  mounted  in  a  naturally  aspirated  ra^ation  shield. 
Lockhart  (1989A)  notes  that  it  is  critically  important  to  maintain  comparable 
aspiration  rates  over  sensors  used  for  temperature  difference  measurements. 

Jacobs  and  McNaughton  (1994)  find  that  a  thin  coat  of  opti^  white  paint  can 
significantly  reduce  solar-radiation  absorption  without  causing  adverse  effects  on 
sensor  response. 

A.3.3  Thermal  Contact.  Probe  wires  and  supports  provide  heat  sotmces  or  smks 
for  a  temperature  sensor.  These  thermometer  components,  like  the  probe  itself, 
absorb  and  emit  radiation.  Heat  transfer  occurs  anytime  there  is  a  temperature 
difference  between  the  probe  and  the  components  to  which  it  is  connected.  Care¬ 
ful  probe  design  and  an  installation  that  thermally  isolates  the  probe  as  much  as 
possible  can  minimize  these  effects. 

A.3.4  Foreign  Subfitanrft  Aeeumulation.  The  accumulation  of  forei^  substances 
or  the  formation  of  oxide  surfaces  on  the  probe  can  change  the  probe’s  response 
time  and  radiative  characteristics.  Probes  should  be  periodically  inspected  and 
cleaned  or  replaced  to  minimize  these  effects. 

A.3.5  Field  Checks.  The  sensor  element  performance  can  be  checked  in  the  field 
through  use  of  one  or  more  isothermal  dewars.  A  dew^  is  characterized  by  a 
well-insulated  core  where  temperature  gradients  are  minimized.  Temperature 
probes  inserted  into  the  dewar  will  quickly  come  into  equihbrium  with  the  interior 
of  the  dewar  and  produce  a  stable  temperature  reading.  This  procedure  is  particu¬ 
larly  useful  for  cb^cking  the  match  of  probes  used  for  temperature  differential  (AT) 
measurements  where  interprobe  functional  precision  is  more  important  than 
absolute  accuracy.  For  this  purpose,  Lockhart  (1989A)  recommends  the  use  of 
several  dewars,  one  chilled  by  ice,  one  kept  to  a  temperature  near  40  C,  and  a 
third  maintained  near  ambient  conditions. 


Exposure  to  electromagnetic  radiation  can  be  minimized  by  isolating  the 
thermometer  unit  as  much  as  possible  firom  other  equipment  and  by  conditioning 
thft  electrical  power.  Long  parallel  cable  runs  should  be  avoided,  particularly  near 
power  cables.  Grounding  should  be  stable  and  common  between  the  sensor  and 
data  recorder. 

A.3.6.  Other  Envjrnnmpnt.al  Stresses.  Other  environmental  stresses  that  impose 
noise  or  variable  biases  on  ambient-air  temperature  measurements  include 
thermal  or  mechanical  shock,  electromagnetic  radiation,  and  self-heating. 
Resistance  temperature  detectors  (RTDs)  include  a  platinum  alloy  film  or  wire 
encased  for  protection  in  a  glass  or  ceramic  matrix.  The  coefficient  of  thermal 
expansion  for  the  metal  element  and  its  matrix  differ,  and  thermal  shock  causes 
stress  on  the  element  that  may  alter  its  cahbration.  All  temperature  probes, 
particularly  quartz  crystal,  are  susceptible  to  mechanical  shock,  which  alters 
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caHbration.  The  RTDs,  thermocouples,  and  thermistors  are  subject  to  electeo- 
magnetic  contamination,  groimd  faults,  or  differences  in  ground  between  the  probe 
and  the  data  recorder.  Noise  from  amplifiers  and  other  system  components  md 
60-Hz  line  noise  may  also  degrade  the  measurement.  Thermistors  and,  to  a  lesser 
extent,  RTDs  cannot  completely  dissipate  heat  generated  by  passage  of  cuirent 
through  the  sensor  elements.  This  excess  heat  produces  a  measurement  bias. 
These  effects  are  minimized  through  use  of  end-to-end  signal  checks,  by  inter¬ 
comparing  probes  immersed  in  an  isothermal  dewar,  and  by  minimizing  exposure 
to  environmental  stresses. 

End-to-end  field-signal  checks  should  detect  system  noise  and  bias  not 
associated  with  the  sensing  element.  These  checks  include  replacing  the  sensor 
vdth  a  precision  decade  resistor  and  measuring  the  output  at  the  data  recorder. 
Biases  can  then  be  zeroed  and  noise  can  be  traced.  Several  sequential  signal 
checks  may  be  needed  to  find  low-frequency  drift  in  the  system. 

A.4  Hygrometers 

Hygrometer  performance  degradation  is  primarily  due  to  contamination  of 
the  sensor  or  control  surface.  Frequent  cleaning,  recahbration,  or  sensor 
replacement  may  be  needed  to  ensure  data  quahty. 

A.4.1  Lithium  Chloride  Dew  Cell.  The  dew  cell  is  a  saturated  salt  dew  point 
instrument  with  a  sensor  consisting  of  a  fabric  bobbin  coated  mth  lithium  chloride 
(LiCl).  An  alternating  current  is  passed  through  the  salt  coating  to  produw 
resistive  heating.  This  heating  causes  evaporation  fimm  the  hygroscopic  LiCl 
solution  with  the  rate  of  evaporation  controlled  by  vapor  pressime.  As  the  bobbin 
loses  water,  resistance  of  the  LiCl  coating  decreases.  The  bobbin  cycles  through 
heating  and  cooling  untfi  equilibrium  with  the  atmospheric  vapor  pressure  is 
reached.  The  eqmlibrium  temperatme  is  used  as  a  measure  of  ambient  dew  point, 
with  an  accuracy  near  ±1  °C  over  a  temperature  span  from  -12  to  +38  °C.  The 
instrument  operates  within  a  range  of  11  to  95  percent  relative  humidity,  but 
performance  declines  markedly  at  humidities  below  20  percent.  Response  is  slow 
(dew  cell  response  time  is  typically  several  minutes)  because  of  bobbin  thermal 
mass  and  the  heating/cooling  cycles;  an  averaging  period  of  at  least  15  minutes  is 
recommended.  The  sensor  shoiild  also  be  operated  within  an  aspirated  radiation 
shield.  Exposiire  to  water  or  excessive  contamination  will  ruin  the  bobbin  coating. 
Once  set  into  operation,  the  instrument  should  be  run  continuously  to  prevent 
moisture  accumulation  from  disturbing  the  LdCl  coating.  Recoating  is  a  simple 
procedure  easily  accomplished  in  the  field  and  should  be  done  weekly  to  minimize 
contamination  effects. 

A.4.2  nbillpd  Mirror  Hygrometer.  Optical  condensation  hygrometers  measime 
dew  point  by  thermoelectric  (Peltier)  cooling  of  a  polished  mirror  smfface  until 
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condensation  occurs.  Dew  formation  on  the  mirror  interrupts  a  light  beam  to  shut 
off  the  cooling  process.  A  thermometer  embedded  in  the  mirror  surface  measures 
the  dew  point.  These  instruments  are  capable  of  operation  over  greater  tempera¬ 
ture  and  humidity  ranges  than  most  other  hygrometers  designed  for  field  use. 

The  response  of  an  optical  condensation  hygrometer  is  limited  by  the  finite 
time  required  to  condense  or  evaporate  a  fiilm  on  the  controlling  mirror  surfa^. 
Accuracies  of  ±0.5  °C  or  better  are  possible.  Operator  adjustments  for  "gain," 
"compensation,"  and  "thickness"  affect  instrument  performance.  Gain  controls  the 
system-dynamic  response,  while  compensation  is  a  phase  lead  to  the  amplifier 
circuit  to  compensate  for  the  phase  lag  of  the  Peltier  cooler.  This  compensation 
dampens  oscillations  to  allow  a  higher  gain  setting  and  to  improve  dynamic  perfor¬ 
mance.  Thickness  controls  the  depth  of  the  film  allowed  on  the  mmor  before 
cooling  is  interrupted.  The  thickness  adjustment  offers  a  compromise  between 
fast  dynamic  response  and  insensitivity  to  mirror  contamination.  Skill  and 
patience  are  needed  to  optimize  these  adjustments,  although  microcontroller 
technology  has  recently  been  used  to  automate  the  adjustment  procedure. 

The  cooling  rate  of  the  Peltier  cooler  in  an  optical  condensation  hygrometer  is 
rapid,  but  settling  time  is  inversely  proportional  to  vapor  pressure.  For  very  low 
vapor  pressures  (usually  associated  with  temperatures  below  -18  °C),  the  small 
amount  of  water  vapor  present  may  require  several  tens  of  minutes  to  form  the 
required  condensate  thickness  on  the  controlling  mirror.  In  spite  of  the  increased 
time  constant,  the  optical  condensation  hygrometer  produces  usefifi  data  at  cold 
temperatures  where  the  response  of  most  other  humidity  sensors  is  poor.  An 
additional  complication  at  temperatures  between  0  and  -30  °C  is  that  either  dew 
or  fi*ost  may  form  on  the  mirror  surface.  The  possibility  of  accumulating  either 
dew  or  fi-ost  creates  a  measurement  ambiguity  unless  fi*eezing  is  artificially 
initiated  on  the  mirror. 

The  mirror  surface  of  the  optical  condensation  hygrometer  is  relatively  insensi¬ 
tive  to  contamination,  although  cleaning  every  2  or  3  days  is  recommended  in  a 
very  dvisty  environment.  Mirror  self-cleaning  procedures  are  generally  ineffective 
in  high  dust  concentrations.  Dirty  mirror  surfaces  can  caxise  the  equipment  to 
oscillate  about  the  condensation  point. 

A.4.3  Adsorption  Hygrometers.  An  adsorption  hygrometer’s  sensor  electrically 
measures  the  percent  relative  humidity  based  on  changes  in  conductivity, 
capacitance,  or  resistance  because  of  adsorption  or  desorption  of  moisture  as 
humidity  changes.  These  hygrometers  typically  operate  over  a  humidity  range  of 
10  to  98  percent.  Soluble  salts,  halogens,  or  any  substance  that  alters  surface 
resistivity  can  decalibrate  or  destroy  this  type  of  sensor.  Consequently,  these 
instruments  should  not  be  used  in  highly  contaminated  envko^ente.  If 
temperature-compensated  and  mounted  in  an  aspirated  radiation  shield. 
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adsorption  hygrometers  can  produce  accuracies  of  ±2  to  5  percent  oyer  a^hmted 
range  of  temperatures.  Major  accuracy  problems  are  likely  at  cold  (-20  C) 
temperatures,  and  significant  calibration  drift  and  sensor  failures  c^  o^  m  a 
marine  environment.  Some  instruments  exhibit  noticeable  hysteresis  effects, 
particularly  following  long  exposure  at  saturation. 

Modem  relative  humidity  (RH)  sensors  are  generally  thin-film  capacitors 
designed  to  change  capacitance  as  a  function  of  moisture  adsorbed  onto  the 
surface  The  primary  source  of  error  is  surface  contamination  which  appears  as  a 
bias  in  the  reading.  Shielding  is  used  to  protect  the  sensor,  but  airborne  contami¬ 
nants  eventuaUy  accumulate  on  exposed  surfaces.  Errors  can  appe^  as  a  slow 
drift  of  cafibration  or  as  an  abrupt  shift  in  the  humidity  v^ue.  Field  checks 
against  an  uncontaminated  instrument  can  be  used  to  verify  performance. 
Immersion  of  the  sensor  into  a  capsule  with  controlled  humidity  is  a  more  ^e 
but  effective  check.  Some  humidity  sensor  manufacturers  provide 
capsules  for  this  purpose  with  humidity  controlled  by  saturated  salt  solutions, 
typically  ranging  firom  20  to  80  percent  RH. 

A.4.4  J.ymftn-n  Hygrometers.  The  Lyman-a  hygrometer  offers  a  fast-response 
fine-scale  water-vapor  fluctuation  measurement  capability.  Instrument  operation 
is  based  on  Beer’s  law  absorption  of  light  at  the  Lyman-a  (121.5-nm)  wavelength 
by  the  hydrogen  constituent  of  water  molecules.  The  sensor  consists  of  a  glow 
discharge  tube  and  a  detector  tube  connected  to  an  electrometer  circmt  with  the 
detector  set  behind  a  magnesium  fluoride  window.  A  gap  of  0.5  to  2  cm  between 
the  source  and  detector  offers  fine  spatial  resolution  at  the  expense  of  signmcant 
flow  distortion  effects.  The  window  is  sensitive  to  etching  and  contamination, 
mgldng'  the  Lyman-tt  subject  to  rapid  cahbration  drift.  Consequently,  the 
Lyman-a  hygrometer  is  often  used  in  conjunction  with  a  slow-respo^e  reference 
instrument  such  as  a  chilled  mirror  hygrometer.  Buck  (1985)  describes  principles 
of  a  Lyman-a  hygrometer,  and  Lind  and  Shaw  (1991)  describe  techmques  for 
adjusting  Lyman-a  cahbration  for  airborne  measurements. 

A.4.5  Tnfrflrpd  Hygrometers.  Infrared  hygrometers  use  closely  spaced 
colunmated  beams  of  infrared  light  and  a  differential  absorption  technique  to 
measure  absolute  humidity.  These  instruments  have  a  transmitter  which  sends 
out  two  parallel  beams  closely  spaced  in  frequency,  but  with  different  water  vapor 
absorption  characteristics.  The  ratio  of  the  energy  received  from  the  absorbing 
and  nonabsorbing  beams  is  proportional  to  the  absolute  humidity  in  the  sampled 
volmne.  With  sealed-optics  and  nonreactive-sensor  plates,  optical-humidity 
sensors  are  relatively  insensitive  to  contamination,  are  capable  of  operation  at 
extreme  temperatures  without  degradation,  and  pro^dde  lugh  data  rates.  Internal 
calibration  checks  minimize  calibration  drift.  Variations  in  chopper-motor  speed 
and  instrument  casing  temperature  can  cause  erroneous  readings. 
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A.4.6  Psvchrometer.  The  psychrometer  consists  of  two  matched  thermometers 
mounted  together  in  a  holder.  One  thermometer  bulb  is  exposed  to  ambient-air 
conditions  (the  dry  bulb)  and  the  other  is  covered  by  a  wet-muslin  sock  (the  wet 
bulb).  When  the  thermometers  are  aspirated,  evaporative  cooling  of  the  wet  bulb 
lowers  its  temperature.  The  difference  between  dry  and  wet  bulb  temperatmes  is 
the  wet-bulb  depression.  Standard  meteorological  tables  relate  dry  bulb  tempera¬ 
ture  and  wet  bulb  depression  to  relative  humidity  and  dew  point.  Some  psychro- 
meters  have  been  automated  but  most  are  used  manually.  Psychrometers  c^ 
provide  very  accurate  readings  for  temperatures  above  0  °C  and  high  humidities. 
At  low  humidities,  a  significant  amount  of  energy  must  be  dissipated  by 
evaporative  cooling  for  the  wet  bulb  to  reach  equibbrium.  At  this  point,  the 
aspiration  rate  becomes  a  crucial  bmitation  on  accuracy.  At  temperatures  below  ^ 

0  °C,  the  wet  bulb  can  be  either  liquid  or  fi*ozen,  each  with  unique  thermodynamic 
equilibrium  values.  Accuracy  is  a  fimction  of  the  care  with  which  the  wet-bulb 
reading  is  taken.  Factors  limiting  accuraQr  include  cleanliness  of  the  muslin  sock, 
water  purity,  and  effects  of  water  temperatme  on  wet-bulb  readings.  Measure¬ 
ment  precision  of  a  liqxaid-in-glass  psychrometer  is  limited  by  the  abifity  of  the 
observer  to  resolve  the  thermometer’s  temperature  scales;  scale  readings  finer 
than  0.1  °C  are  probably  not  meaningful.  The  dynamic  response  of  a  psychrom¬ 
eter  is  usually  slow,  although  Tsxikamoto  (1986)  describes  a  fast  fine-wire 
psychrometer  used  for  water  vapor  flux  measurements. 

A.5  Actinometers 

Actinometers  should  be  mounted  level  (or  at  a  specific  chosen  angle)  on  a 
rigid  surface  fi*ee  of  obstructions  above  the  plane  of  the  sensing  element. 

A^5.1  Foreign  Rnhstarirfi  AcciiTnulation  on  Dome.  Exposed  sensor  surfaces  or 
domes  are  susceptible  to  foreign  substance  accumulation.  These  surfaces  require 
periodic  inspection  and  cleaning. 

A.5.2  Calibration  Drift.  Exposure  of  the  instruments  to  wide  temperature 
changes  or  thermal  gradients  may  cause  the  instrument  to  drift  off  calibration. 

A.5.3  WnTniditv.  Moisture  trapped  inside  the  dome  covering  the  sensor  causes 
erroneous  readings  and  is  minimized  through  the  use  of  a  desiccant.  The 
desiccant  should  be  inspected  and  changed  as  necessary  to  minimize  moisture 
accumiilation. 

A.5.4  Exnosure  Blockage.  Care  shoxild  be  taken  to  avoid  inadvertent  shadowing 
of  solar  radiation  sensors.  The  instrument  site  should  be  selected  to  ensure  that 
trees,  bmldings,  or  other  structure  do  not  cast  a  shadow  on  the  instrument  during 
any  part  of  the  day  when  data  are  required.  Reflections  fi-om  nearby  objects  such 
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as  windows  or  bright  sheet  metal  should  be  avoided.  It  is  desirable  to  <^ose  a 
site  where  all  obstades  subtend  an  angle  of  less  than  10  above  the  horizon. 


A-6  Barometer/Altimeter 


Some  modem  pressure  sensors  are  hand-held  devices  that  serve  as 
altimeters  as  well  as  pressure  sensors.  These  instruments  are  desired  for  fie  d 
use  and  are  usually  temperature  compensated  to  provide  accurate  altimetry. 

Older  style  barometers  of  the  mercury-in-glass  and  aneroid  type  do  not  travel  well 
and  should  be  kept  in  stable  environments  where  they  are  not  subject  to  abrupt 
temperature  or  pressure  changes. 


Exposure  is  a  major  consideration  for  barometers,  especially  those  located 
inside  buildings.  Indoor  pressure  changes  occur  as  doors  open  or  dose  and  as  ^ 
conditioning  units  cycle  on  or  off.  Exposure  to  outside  pressure  can  be  established 
through  an  orifice,  but  the  opening  can  become  blocked  by  debns  or  damage  to  the 
connecting  tube.  Building  wake  effects  may  influence  pressure  readmgs  durmg 
high  winds.  Mechanical  shock  can  decalibrate  an  instrument. 


A.7  Aerosol  Samnlers 

Unlike  atmospheric  gases,  it  cannot  be  assumed  that  aerosols  ^spended  in 

the  air  will  be  well  mixed.  Aerosol  behavior  is  governed  by  inertial  effects, 

settling  velodties,  electrostatic  effects,  and  shape  factors.  The  aerosol  conteiff  of 
the  air  is  measured  by  inertial  impaction,  dosimetry,  and  by  nephelometry.  Each 
of  these  measurement  techniques  is  affected  by  the  characteristics  of  the  aerosol  to 
be  measured.  Although  calibration  of  aerosol  samplers  is  relatively  straight¬ 
forward  and  subject  to  few  fault  conditions,  obtaining  a  sample  that  is  truly 
representative  of  the  ambient-aerosol  distribution  can  be  difficult.  Aerosols  are 
found  on  scales  ranging  fi-om  less  than  one  micron  to  over  100  microns  m  diameter 
with  behaviors  that  are  dependent  on  their  size,  shape,  and  chenucal  makeup. 
Unrepresentative  aerosol  sampling  typically  occurs  for  the  following  reasons: 

1.  samplers  located  where  they  are  inadequately  exposed  to  the  aerosol  or 
are  blocked  by  other  objects; 

2.  mismatch  between  wind  speed  and  sampler  flow  velocity  or  misalign¬ 
ment  of  sampler  in  the  flow  (nonisokinetic  effects); 


3.  lack  of  pre-  and  post- test  baselines; 

4.  sampling  time  so  short  that  an  insiiffident  sample  is  obtained  for 
analysis,  or  so  long  that  filters,  optics,  or  inlets  become  overloaded  or  blocked; 
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5.  particle  shape  or  electrostatic  effects  skew  the  sample;  and 

6.  aerosol  cloud  is  too  thin  or  too  dense  for  the  chosen  sampler. 

A.7.1  Inertial  Imnactors.  Inertial  impactors  are  used  to  obtain  mass  measure¬ 
ments  of  aerosol  content,  usually  with  the  total  mass  separated  into  size  ranges. 
The  performance  of  inertial  impactors  is  dependent  on  inertial  particle-size 
separation  as  the  aerosol  doud  passes  over  various  flow  obstructions.  A  particle’s 
behavior  is  assumed  to  vary  with  its  equivalent  diameter,  adhering  to  collector 
plates  by  electrostatic  and  Van  der  Waals  forces.  Larger  and  denser  particles 
t3q)ically  are  the  first  to  contact  the  impactor,  although  impaction  and  separation 
is  a  stochastic  process  with  an  overlap  in  size  ranges.  Collection  efficiency  vanes 
with  particle  size,  shape,  density,  and  electrostatic  charge.  Particles  that  have 
facets  or  points  tend  to  accumulate  charge,  while  round  particles  are  less  likely  to 
retain  electric  charges.  Atmospheric  moisture  content  also  affects  collection 
efficiency.  Inertial  impactors  should  be  operated  for  a  sufficient  length  of  time  to 
accumulate  a  large  enough  mass  for  measurement.  Measurement  accuracy  is 
dependent  on  how  accurately  the  accumulated  particle  masses  can  be  weighed. 

A.7.2  noRimfttftrfi.  Dosimeters  include  filters  that  collect  particles  by  diffusion, 
inertial  separation,  and  direct  interception;  electrical  charges  on  the  particles  and 
filters  generally  increases  collection  efficiency.  If  sampling  is  continued  too  long, 
filter  collection  efficiency  decreases  as  aerosols  b3q)ass  the  clogged  filter.  Dosage  is 
calculated  as  a  function  of  the  mass  of  material  collected,  divided  by  volumetric 
flowrate.  (NOTE:  The  effects  of  measurement  site  pressures  must  be  considered 
if  volumetric  flow  rates  taken  at  different  sites  are  to  be  compared.)  Dosimeters 
obtain  their  most  representative  samples  when  the  flow  rate  is  matched  with  the 
speed  at  which  the  wind  propels  aerosols  towards  the  filter.  Representativeness  is 
therefore  a  function  of  inlet  geometry,  orientation  with  respect  to  the  wind,  and 
the  degree  to  which  the  flow  rate  is  matched  to  wind  speed  (isokinetic  sampling). 

It  is  difficult  to  obtain  representative  dosimeter  measurements  for  particle  sizes 
larger  than  25  microns. 

A.7.3  Nenhelometers.  Nephelometers  operate  by  nonextractive  optical 
techniques,  defining  aerosol  density  using  light  scattering  or  laser-Doppler 
velocimetiy.  Nephelometer  calibration  is  performed  using  standard  polystyrene 
latex  beads  available  from  National  Institute  of  Standards  and  Technology  (NIST) 
or  commercial  sources.  Light-scattering  systems  have  measurement  uncertainties 
because  of  nonspherical  particle  shapes  and  aerosol  refi*active-index  variations. 
Assuming  constant  particle  shape  and  density,  laser  velodmetry  uses  particle- 
velodty  measurements  in  an  accelerating  flow  field  to  define  size.  Nephelometry 
produces  unrepresentative  readings  if  the  aerosol  is  so  dense  that  multiple 
scattering  of  the  light  occurs,  if  particle  collisions  alter  the  size  distribution  during 
the  sampling  period,  or  if  the  aerosol  accumulates  on  the  system  optics. 
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A.8  Scintillometers 


A  scintillometer  is  a  ground-based  remote-sensing  insti^ent  desired  to 
measure  crosswind  components  and  optical  turbulence  intensity  along  a  line- 
of-sight  path  established  between  a  transmitter  and  a  dow^ange  receiver  (see 
Eange  Commanders  Council  document  381-92,  Meteorological  Measurements 
Guide).  Typical  scintillometer  operating  range  is  0  to  20  m/s  for  crosswinds  and 
10'^®  to  10'^^  for  the  refractive  index  structure  parameter  (C^^).  The  pathlength 
between  transmitter  and  receiver  can  range  from  a  few  hundred  meters  to  several 
kilometers,  depending  on  the  along-path  turbulence  intensity.  Excessive  turbu¬ 
lence  causes  signal  saturation,  where  the  relationship  between  the  log  ampHtude 
of  the  variance  of  signal  intensity  and  C„^  is  no  longer  linear.  Newer  scintil¬ 
lometer  models  can  compensate  for  some  intrusion  into  the  nonlinear  region,  but 
saturation  should  he  avoided  to  preserve  data  quality.  The  theory  and  perfor¬ 
mance  of  a  recent  model  crosswind  profiling  scintillometer  is  described  by  Ochs  et 
al.  (1992). 

Because  scintillometers  sense  changes  in  the  intensity  and  position  of 
optical  turbulence  as  the  basis  for  measurement,  they  are  sensitive  to  platform 
vibrations  and  must  be  mounted  on  sohd,  vibration-free  platforms.  Transmitter 
and  receiver  optics  must  also  be  protected  from  direct  mtrusion  of  sunhght. 
Optical-path  positions  must  be  carefully  chosen  or  shaded  to  avoid  direct  sun,  or 
damage  may  occur.  Scintillometer  transmitters  and  receivers  can  he  set  up  to 
operate  inside  buildings  through  windows  as  long  as  excessive  glare  or  reflected 
light  does  not  enter  the  optics.  Scintillometer  signals  are  modulated  to  minimize 
interference  from  backgroxmd  movement,  but  excessive  movement  in  the  receiver's 
field  of  view  can  induce  noise.  For  example,  the  optical  path  should  not  be 
established  where  road  traffic  is  moving  behind  the  transmitter.  Objects  moving 
across  the  optical  path  between  the  transmitter  and  receiver  causes  signal  loss. 

A.9  Sonic  Anemometer/Thermometers 

A  sonic  anemometer/thermometer  consists  of  a  transducer  array  containing 
paired  sets  of  acoustic  transmitters  and  receivers,  a  system  clock,  ^d  micro¬ 
processor  circuitry  for  measuring  intervals  of  time  between  transmission  and 
reception  of  sound  pulses.  Wind  velocity  and  the  speed  of  sound  can  be  calculated 
from  measurements  of  the  transit  time  of  an  acoustic  wavefront  along  the  known 
path  between  pairs  of  transmitter/receivers.  Measurements  are  made  along  one, 
two,  or  three  axes  to  define  components  of  the  three-dimensional  windfield. 
Measurement  accuracy  is  a  function  of  how  precisely  time  is  measured  and  how 
precisely  the  acoustic  pathlength  is  known.  Factors  adversely  affecting  measure¬ 
ment  accuracy  include  changes  in  acoustic  pathlength  caused  by  mechamcal  or 
thermal  stress,  distortion  or  attenuation  of  the  acoustic  wavefront,  or  changes  in 
the  transducers  or  electronic  circuitry  (decalibration).  Flow  distortion  and 
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blockage  by  the  transducers  or  the  supporting  array  can  adversely  affect  the 
velocity  component  measiorement  representativeness. 


Sonic  anemometer/thermometer  fault  conditions  occur  when  tlie  awustic 
wavefront  is  distorted  or  blocked,  when  the  t^durers  or  elecfromcs  M  or 
become  unstable,  or  when  calibration  is  lost.  The  ultr^omc  pulse  emtted  by  the 
transmitter  is  sensitive  to  environmental  interference  (ne^y  explosion  or  jet  am- 
craft),  to  acoustic  reflections  from  objects  placed  too  dose  to  the  acoustic  may,  or 
to  attenuation  or  blockage  by  predpitation,  dust  p^des,  or  inse^  on  the  fr^- 
ducers.  Nonreceipt  of  an  acoustic  wavefront  usually  causes  a  no-^  , 

Noise  or  intermittent  transducer  failures  appear  as  noise  spikes.  Themal  md 
shock  may  cause  decalibration.  (Cahbration  can  be  verified  m  the 
field  through  use  of  a  zero-wind  chamber.) 


A.10  r.nnfitant  Temperature  Anemometers 

Constant  temperature  anemometer  (CTA)  systems  consist  of  ^plifiers,  a 
data-acquisition  system,  and  a  hot-wire  or  hot-film  element  mounted  in  a 
Wheatstone  bridge  that  is  balanced  to  maintain  the  sensor  at  a  conteol  tempera- 
ture  well  above  the  ambient  temperature.  As  wind  blowing  across 
the  heated  element,  the  bridge  network  responds  by  sending  ciment 
element  until  balance  is  restored.  System  electronics  sense  and  record  toe 
resulting  voltage  changes  across  toe  bridge.  The  heat  loss  eiyenenced  by  t  e 
element  (parameterized  by  the  Nusselt  number  Nu)  is  related  to  toe  flow  velocity 
(parameterized  by  Reynolds  number  Re)  across  toe  element  by 


Nu  =  A  +  B  Re^ 


(A-1) 


where  A  and  B  are  empirical  constants  and  p  is  a  velocity-dependent  exponent 
varying  firom  0.45  to  0.51. 


Misorientation  and  contamination  are  the  major  sources  of  error  with  a 
CTA.  The  constant-temperatiu-e  element  is  most  sensitive  to  the  flow  component 
normal  to  it.  Except  for  omnidirectional  sensors,  toe  element  should  be  care^ly 
aligned  into  toe  flow  to  be  measured.  Caps  and  baffles  ^e  o^n  used  to  protect 
the  sensor  fi-om  moisture  and  contamination.  This  housing  distorts  toe  flow,  and 
eddies  shed  from  probe  support  structures  can  cause  unrepresentafive  readings. 
Accumulation  of  rime  ice,  snow,  or  dust  on  toe  sensor  causes  decahbration,  ^d 
accumulations  on  the  sensor  housing  increase  flow  blockage.  When  operated  m 
high-speed  flows,  a  sensor  may  vibrate.  Vibration  causes  a  sixain  on  toe  elemen  s 
which  leads  to  changes  in  wire  resistance.  The  use  of  low-resistance  wes  or  films 
and  low-voltage  signals  also  makes  these  instruments  sensitive  to  electromagnetic- 
field  changes,  electronic  noise,  power-supply  noise,  and  zeroing  drift. 
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Because  convective  cooUng  on  the  surface  of  a  CTA  sensor  is  a  molecular 
heat-exchange  process,  the  relationship  between  cooling  and  wind  speed  should  be 
density  compensated  (compensated  for  changes  in  both  temperature  and  pressure). 
Therefore,  a  (DTA  system  calibrated  at  standard  temperature  and  pressure  (STP) 
would  not  be  in  calibration  for  operation  at  conditions  that  differ  from  STP 
Because  of  the  nonlinear  relationship  between  velocity  and  coohng,  calibrations 
should  mclude  several  points  spanning  the  intended  operating  range  of  the 
instrument.  Practices  for  using  CTA  systems  for  ground-based  and  aircraft 
measurements  are  described  by  Lenschow  (1986). 

A.11  Field  Mill 

electromechanical  devices  designed  to  measure  the  electric 
field  over  flat  surfaces  where  the  fair  weather  potential  gradient  is  uniform.  The 
insists  of  a  rotating,  grounded,  fanlike  plate  that  chops  the  electric  field 
Imes  by  alternately  covering  and  imcovering  plates  called  stators.  The  stators 
charge  and  discharge  in  proportion  to  the  amphtude  of  the  electric  field,  and  the 
de^^^  charge  is  converted  to  voltage  for  output  to  a  recording 

™  characteristics,  contamination,  and  mechanical  feilure  are 
the  t^ee  major  sources  of  error  for  field  mills.  If  a  field  mill  is  located  near  tall 
objects,  the  electric-field  amplitude  will  be  altered.  A  mill  surrounded  by  a  25-foot 
ra^us  of  gravel,  a  50-foot  radius  of  cleared  vegetation,  and  with  no  obstructions 
su^nding  an  angle  of  greater  than  18°  above  the  horizon  is  likely  to  have  a 
s}™ciently  luuform  fair  weather  potential  gradient.  Other  site-dependent 
ch^actenstics  mclude  the  proximity  of  sources  of  charged  particles  such  as 
vemcular  exhai^t,  nearby  power  lines,  or  blowing  sand  or  dust.  Roughening  of 
p  ate  surfaces  that  oreurs  with  corrosion  or  contamination  causes  premature  plate 
discharge,  plate  msulation,  or  physical  contact  between  the  rotor  and  stator.  Field 
mms  are  not  very  sensitive  to  changes  in  rotor  speed  and  will  not  fimction  if 
rotation  ceases. 

Field  m^s  used  in  a  network  are  usually  mounted  and  calibrated  with  the 
same  geometric  and  electrical  characteristics.  Mounting  the  miU  in  an  inverted 
w^guration  ^mizes  the  effect  of  rain  or  particulate  deposition  on  the  plates 
men  operated  m  a  network,  a  comparison  of  the  outputs  of  adjacent  mills  can  be 
®  ofv^  performance.  Systematic  differences  between  nearby 

^Is  ^at  persist  for  an  hour  or  more  are  suspect  and  should  be  investigated 
After  heay  storm  activity,  mill  calibration  should  be  performed  to  verify  that  the 
sensor  and  support  electronics  were  not  damaged  by  lightning. 
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APPENDIX  B 

INSTRUMENT  QUALITY  CONTROL  CHECKS 
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FLAGGING  CRITERIA 


This  appendix  contains  flagging  criteria  with  nominal  numerical  values  that 
pertain  to  a  mid-latitude  site  with  open  exposiire  in  flat  terrain  and  nonextreme 
weather  conditions.  Some  criteria  wiU  work  better  at  certain  sites  than  at  others. 
Niunerical  values  should  be  adjusted  for  site-specific  conditions. 

B.l  R/>t.ating  Anemometer  QC  Checks 


FLAGGING  CRITERIA  FOR  ROTATING  ANEMOMETERS 

Scalar  wind  speed  <0.1  or  greater  than  19.9  m/s  (or  similar  user-defined 
Tr.a-dTmim-TniniTnTiTn  criteria  derived  fi*om  equipment  operating  range)* 

Wind  speed  >2  m/s  and  increases  or  decreases  by  a  factor  of  2  or  greater 
between  wind  readings  (inadequate  sensor  exposure). 

Wind  speed  makes  a  step  increase  fi'om  calm  to  >1  m/s.  (Worn  bearings 
or  ice-bound  sensor  is  likely.) 

Wind  speed  remains  calm  or  at  a  constant  value  for  successive  averaging 
periods  up  to  1  hoiur  (sensor  stuck  in  position). 

Point  on  wind-speed  profile  deviates  fi'om  profile  mean  by  >0.5  m/s 
(possible  bearing  drag). 

Battery  voltage  drops  below  system  reference  voltage. 

Difference  between  TnarimiiTn  and  minimum  wind  speeds  exceeds  the 
average  wind  speed  by  a  factor  of  2.  (Average  possibly  contaminated  by 
noise  spikes  or  signal  drop-out.) 


♦Because  of  the  nonlinear  response  of  rotating  anemometers  at  low  speeds,  the 
calibration  curve  for  some  instruments  may  cause  a  small  negative  wind 
reading  (-0.1,  for  example)  during  calm  winds.  Flagging  criteria  may  require 
adjustment  to  accommodate  this  condition. _ _ _ _ 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 
(7) 
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B.2  Wind  Vane  QC  Checks 


FLAGGING  CRITERIA  FOR  WIND  VANES 

(1)  Wind  direction  remains  constant  for  1  hour.  (Vane  is  stuck  in  position.) 

(2)  Wind  speed  is  >2  m/s  and  wind-direction  standard  deviation  is  less  than 
1°  or  greater  than  60°.  (Vane  is  stuck  in  position  or  inadequately 
exposed.) 

(3)  Wind  speed  is  >2  m/s  and  wind  direction  changes  by  more  than  45° 
between  successive  readings  (inadequate  vane  exposure). 

(4)  Wind  speed  is  >2  m/s  and  point  on  wind-direction  profile  deviates  fi-om 
profile  mean  by  more  than  45°. 

(5)  The  modal  (most  fi-equent)  wind  direction  reading  does  not  change 
between  successive  averages  (dead  spot  on  potentiometer). 

(6)  Battery  voltage  drops  below  system  reference  value. 
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B.3  Thermometer  QC  Checks 


FLAGGING  CRITERIA  FOR  THERMOMETERS 

(1)  Temperature  remains  invariant  or  changes  less  than  0.1  °C  within  a  15-60 
minute  averaging  period. 

(2)  Temperatxire  readings  exceed  dimatological  extremes  (minima  or  maxima) 
for  the  measurement  site  (inadequate  sensor  exposure). 

(3)  Temperature  change  in  excess  of  10  °C  per  hour  (inadequate  sensor 
exposure). 

(4)  Data  point  deviates  from  temperature  profile  by  >0.5  ®C  (for  mechanically 
aspirated  AT  system)  (aspiration  motor  failure  or  calibration  drift). 

(5)  Temperature  standard  deviation  exceeds  2  °C  for  an  averaging  period  of  1 
horn  or  less  (possible  spikes  or  inadequate  sensor  exposure). 

(6)  Difference  between  temperature  maximum  and  minimum  for  a  15-60 
TniniTnuTn  averaging  period  is  zero  or  greater  than  10  “C  (sensor  failure  or 

spikes). 

(7)  Persistent  autoconvective  lapse  rate  is  present  in  temperature  profile. 
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FLAGGING  CRITERIA  FOR  DEW  POINT  INSTRUMENTS 

(1)  Dgw  point  equnls  or  exceeds  companion  temperature  measurement 
(inadequate  sensor  exposure/no  evaporative  cooling). 

(2)  Dew  point  exceeds  expected  climatological  maximum  or  minimum 
(inadequate  sensor  exposure). 

(3)  Dew  point  change  exceeds  10  “C  within  1  hour  (aspiration  failure  or 
sensor  contamination). 

(4)  Dew  point  standard  deviation  exceeds  1  “C  (spikes  in  data). 

(5)  Difference  between  dew  point  maximum  and  minimum  is  zero  or  greater 
than  10  “C  over  the  averaging  period  (spikes  in  data). 

(6)  Dew  point  oscillation  (control  surface  contaminated). 
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B.5  Hiitnldity  Tnstnimcnt  QC  Checks 

FLAGGING  CRITERIA  FOR  RELATIVE  HUMIDITY  INSTRUMENTS 

(1)  Humidity  reads  ^  or  ^100  percent  (inoperative  sensor  or  no  signal). 

(2)  Humidity  reading  remains  constant  over  a  period  of  2  or  more  hours 
(condensation  or  sensor  or  hysteresis  effects). 

(3)  Humidity  changes  by  30  percent  in  1  hour  or  less  (sensor  contamination) 
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FLAGGING  CRITERIA  FOR  PRESSURE  INSTRUMENTS 


(1)  Pressiire  difference  from  dimatological  mean  station  pressure  by  ±25  mb 
(calibration  error). 

(2)  Pressure  change  exceeds  3  mb  in  one  hour  (inadequate  exposure). 

(3)  Pressure  change  of  0.1  mb  or  less  over  a  24-hour  period  (sensor  failure  or 
inadequate  exposure). 

(4)  Difference  between  nriflYiTniim  and  minimum  pressure  over  a  24-hour 
period  exceeds  15  mb  (inadequate  exposure). 
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B.7  Antinometftr  QC  Checks 


FLAGGING  CRITERIA  FOR  ACTINOMETERS 


(1) 

(2) 

(3) 

(4) 

(5) 


Incoming  solar  radiation  exceeds  solar  constant 
(~1400  W  m'®)  (calibration  drift  or  inadequate  exposure). 

Zero  incoming  radiation  during  daylight  hours  (sensor  failure). 


Incoming  radiation  exceeds  100  W  m'^  at  night  (calibration  drift). 


No  change  in  radiation  rea<^g  over  a  period  of  2  or  more  hours  during 
daylight  hours  (loss  of  sensitivity). 


A  drop  in  radiation  readings  not  associated  with  sunrise  or  sunset  jkotits 
at  approrimately  the  same  time  of  day  for  3  or  more  successive  days 
(nossible  partial  sensor  shadowing). 


B.8  VisibiUtv  Sensor  QC  Checks 


FLAGGING  CRITERIA  FOR  VISIBILITY  SENSORS 

(1)  Unexplained  zero  or  full  scale  readings  for  several  successive  hours 
(instrument  failure). 

(2)  Visibility  reading  includes  intermittent  sharp  spikes  (intermittent  optical 
path  interruption  or  loss  of  signal). 
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FLAGGING  CRITERIA  FOR  AEROSOL  SAMPLERS 

(1)  Particle  distribution  is  monodisperse  (calibration  inadequate). 

(2)  Background  readings  are  unstable  and  include  sharp  spikes  (electronics 
fault). 


B-11 


B.ll  rnngtnnt  Temperature  Anemometers 

FLAGGING  CRITERIA  FOR  CONSTANT  TEMPERATURE  ANEMOMETERS 

(1)  Sudden  "ping"  or  change  in  signal  level  (probe  contamination). 

(2)  Increase  in  noise  level  (electromagnetic  contamination  or  electronics  fault) 

(3)  Suspect  changes  in  signal  characteristics  (straining  of  element  caused  by 
vibration,  change  in  probe  orientation  with  respect  to  wmd,  or  probe 
housing  blockage). 
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B.12  Field  Mills 


FLAGGING  CRITERIA  FOR  FIELD  MILLS 

(1)  Constant  large  readings  during  fair-weather  conditions  (contamination  of 
the  mill  or  a  support  electronics  fault). 

(2)  Sensor  data  that  appear  stormlike  during  fair  weather  conditions. 
(Blowing  dust  and  sand  contain  charged  particles  that  can  produce  a 
stormlike  response.) 

NOTE:  The  mill  motor  must  be  operating  to  produce  data.  If  a  motor 
fault  bit  is  available  from  the  mill  support  electronics,  then  the  operating 
software  should  use  the  state  of  this  bit  to  determine  the  validity  of  the 
data.  If  no  fault  bit  is  available,  the  sensor  should  be  checked  daily  for 
motor  operation. 
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